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Abstract 

These lectures describe CP violation and rare decays of K and B mesons and consist of 
ten chapters: i) Grand view of the field including CKM matrix and the unitarity triangle, 
ii) General aspects of the theoretical framework, iii) Particle-antiparticle mixing and CP 
violation, iv) Standard analysis of the unitarity triangle, v) The ratio e'/e including most 
recent developments, vi) Rare decays — > ■k'^vu and n^ui), vii) Express review of 

other rare decays, viii) Express review of CP violation in B decays, ix) A brief look beyond 
the Standard Model including connections between e'/e and CP violating rare K decays, x) 
Final messages. 
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1 Grand View 



1.1 Preface 

CP violation and rare decays of K and B mesons play an important role in the tests of the 
Standard Model and of its extensions. The prime examples, already observed experimentally, 
are — and — mixings, CP violation in Kl — > vrvr and the rare decays B X7, 
Kl — > and ir^vu. In the coming years CP violation in B decays, B^ — B^ mixing 

and rare decays Ki tt^i^u, Ki Tr'^e+e", B — > Xg^J,^!^ , B^g l^l^ and B Xg^^vu 
will hopefully be included in this list. 

These lectures provide a non-technical description of this fascinating field. There is un- 
avoidably an overlap with my Les Houches lectures Q and with the reviews Q and 1^]. On 
the other hand new developments are included and all numerical results updated. 

1.2 Some Facts about the Standard Model 

Throughout these lectures we will dominantly work in the context of the Standard Model with 
three generations of quarks and leptons and the interactions described by the gauge group 
SU{2>)c SU{2)l U{1)y spontaneously broken to SU{3)c U{1)q. There are excellent 
text books on the dynamics of the Standard Model. Let us therfore collect here only those 
ingredients of this model which are fundamental for the subject of weak decays. 

• The strong interactions are mediated by eight gluons Ga, the electroweak interactions 
by W^, Z° and 7. 

• Concerning Electroweak Interactions, the left-handed leptons and quarks are put into 
SU{2)l doublets: 




c 

L \ / L \ / L 




(1.1) 



(1.2) 



with the corresponding right-handed fields transforming as singlets under SU (2)^. The 
primes in (|l.2| ) will be discussed in a moment. 



The charged current processes mediated by are flavour violating with the strength 
of violation given by the gauge coupling g2 and effectively at low energies by the Fermi 
constant 

Gp §2 

and a unitary 3x3 CKM matrix. 
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The CKM matrix || connects the weak eigenstates {d',s',b') and the corresponding 
mass eigenstates d,s,b through 
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(1.4) 



In the leptonic sector the analogous mixing matrix is a unit matrix due to the mass- 
lessness of neutrinos in the Standard Model. 

• The unitarity of the CKM matrix assures the absence of flavour changing neutral current 
transitions at the tree level. This means that the elementary vertices involving neutral 
gauge bosons {Ga, Z^, 7) are flavour conserving. This property is known under the 
name of GIM mechanism [Q. 

• The fact that the V^^'s can a priori be complex numbers allows CP violation in the 
Standard Model |]. 

1.3 CKM Matrix 
1.3.1 General Remarks 

We know from the text books that the CKM matrix can be parametrized by three angles and 
a single complex phase. This phase leading to an imaginary part of the CKM matrix is a 
necessary ingredient to describe CP violation within the framework of the Standard Model. 

Many parametrizations of the CKM matrix have been proposed in the literature. We will 
use two parametrizations in these lectures: the standard parametrization ||^] recommended 
by the Particle Data Group and the Wolfenstein parametrization 



1.3.2 Standard Parametrization 

With Cij = cos Oij and Sij = sinOij = 1,2,3), the standard parametrization is given by: 
^ C12C13 S12C13 sise^'^ \ 

S12C23 - Cl2S23Sl3e"^ C12C23 - Sl2S23Sl3e*'^ S23C13 

C23C13 J 



VcKM 



\ S12S23 - Cl2C23Sl3f 



-S23C12 - Sl2C23Sl3e 



(1.5) 



where 5 is the phase necessary for CP violation. Cjj and Sij can all be chosen to be positive 
and 6 may vary in the range < (5 < 2ii. However, the measurements of CP violation in K 
decays force 5 to be in the range < 5 < vr. 
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From phenomenological applications we know that S13 and ,523 are small numbers: O(10~'^) 
and 0(10^^), respectively. Consequently to an excellent accuracy C13 = C23 = 1 and the four 
independent parameters are given as 

Sl2 = \Vus\, Si3 = \Vub\, S23 = \Vcb\, 5. (1.6) 

The first three can be extracted from tree level decays mediated by the transitions s ^ u, 
b ^ u and b ^ c respectively. The phase 6 can be extracted from CP violating transitions 
or loop processes sensitive to \Vtd\- The latter fact is based on the observation that for 
< 6 < TT, as required by the analysis of CP violation in the K system, there is a one-to-one 
correspondence between 6 and \ Vtd\ given by 

\Vtd\ = Va^ + b2-2abcos6, a = \VM, 6=|KdKfe|- (1.7) 

The main phenomenological advantages of (^^) over other parametrizations proposed in 
the literature are basically these two: 

• si2; si3 and S23 being related in a very simple way to \Vus\, \Vub\ and \Vcb\ respectively, 
can be measured independently in three decays. 

• The CP violating phase is always multiplied by the very small S13. This shows clearly 
the suppression of CP violation independently of the actual size of 5. 

For numerical evaluations the use of the standard parametrization is strongly recom- 
mended. However once the four parameters in ( |1.6P have been determined it is often useful 
to make a change of basic parameters in order to see the structure of the result more trans- 
parently. This brings us to the Wolfenstein parametrization 0] and its generalization given 



in lOl. 



1.3.3 Wolfenstein Parameterization 

The Wolfenstein parametrization is an approximate parametrization of the CKM matrix in 
which each element is expanded as a power series in the small parameter A = \ Vus\ = 0.22, 



/ 1_A! 

^ 2 



V 



A AX^{g-ir]) \ 



-A 1 - ^ AX^ 



V 



+ 0(A^), (1.^ 



AX^{1- Q-i-q) -AX"^ 1 
and the set (|1.6|) is replaced by 

A, A, Q, rj. (1.9) 



3 



Because of the smallness of A and the fact that for each element the expansion parameter 
is actuaUy A^, it is sufficient to keep only the first few terms in this expansion. 

The Wolfenstein parametrization is certainly more transparent than the standard parametriza- 
tion. However, if one requires sufficient level of accuracy, the higher order terms in A have 
to be included in phenomenological applications. This can be done in many ways. The point 
is that since ( |1.8| ) is only an approximation the exact definiton of the parameters in ( |1.9| ) 
is not unique by terms of the neglected order O(A^). This situation is familiar from any 
perturbative expansion, where different definitions of expansion parameters (coupling con- 
stants) are possible. This is also the reason why in different papers in the literature different 



0{X'^) terms in (1.8) can be found. They simply correspond to different definitions of the 
parameters in ( |1.9D . Since the physics does not depend on a particular definition, it is useful 
to make a choice for which the transparency of the original Wolfenstein parametrization is 
not lost. Here we present one way of achieving this. 

1.3.4 Wolfenstein Parametrization beyond LO 

An efficient and systematic way of finding higher order terms in A is to go back to the standard 
parametrization ( pT^ ) and to define the parameters (A, A, g, rj) through [|lO|, |ll| 

si2 = A, S23 = A\\ sne-'^ = A\\q - ir]) (1.10) 

to all orders in A. It follows that 

■coso, r] = smo. (1-11) 



S12S23 S12S23 
(|1.10| ) and ( 1.11 ) represent simply the change of variables from (1^) to (|1.9| ). Making this 



change of variables in the standard parametrization ( |1.5| ) we find the CKM matrix as a 
function of (A, A, g, rf) which satisfies unitarity exactly. Expanding next each element in 
powers of A we recover the matrix in ( |1.S| ) and in addition find explicit corrections of O(A^) 
and higher order terms:. 

= 1 - ^A^ - Ia^ + 0(a6) (1.12) 
Vus = \ + 0{\'), Vub = AX^{Q-ir^) (1.13) 
Vcd = -A + \a^X''[1 - 2{q + i77)] + O(A^) (1.14) 

Vcs = l- - + 4^2) + 0{\^) (1.15) 

Z o 

V^ = A\^ + 0{\''), Vtb = l-\A^X^ + 0{\^) (1.16) 



4 



Vu = AX-" 



{Q + ir}){l-\x^) 



+ 0{X^) (1.17) 



Vts = -AX^ + ^A(l - 2^)A^ - iriAX'^ + 0{X^) (1.18) 

We note that by definition V^h remains unchanged and the corrections to Vus and Vch 
appear only at 0{X'') and 0(A®), respectively. Consequently to an an excellent accuracy we 
have: 

Vus = A, Vcb = AX^, (1.19) 

Vub = AX^{q - irj), Vtd = AX^'il -Q-if)) (1.20) 

with ^ ^ 

Q=q{1-\), f/ = r?(l-y). (1.21) 

The advantage of this generalization of the Wolfenstein parametrization over other gener- 
alizations found in the literature is the absence of relevant corrections to Vus, Vcb and Vub 
and an elegant change in Vtd which allows a simple generalization of the so-called unitarity 
triangle beyond LO. 

Finally let us collect useful approximate analytic expressions for Aj = VidV*^ with i = c,t: 
ImXt = -ImXc = ry^^A^ =| Kfe \ \Vcb\ sin 5 (1.22) 

\2 

ReXc = -A(l - —) (1.23) 

ReXt = -{l-^)A^XHl-Q). (1.24) 

Expressions ( |1.22| ) and ( 1.23| ) represent to an accuracy of 0.2% the exact formulae obtained 
using ( p..5[ ). The expression ( |1.24| ) deviates by at most 2% from the exact formula in the 
full range of parameters considered. For g close to zero this deviation is below 1%. After 
inserting the expressions ( |1.22| )-( [l.24| ) in the exact formulae for quantities of interest, a 
further expansion in A should not be made. 



1.3.5 Unitarity Triangle 

The unitarity of the CKM-matrix implies various relations between its elements. In particular, 
we have 

VudV:b + VM + VM = 0. (1.25) 

Phenomenologically this relation is very interesting as it involves simultaneously the elements 
Vubi Vcb and Vtd which are under extensive discussion at present. 

The relation ( |1.25D can be represented as a "unitarity" triangle in the complex (g, ff) plane. 
The invariance of ( |1.25| ) under any phase-transformations implies that the corresponding 
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triangle is rotated in the {g, fj) plane under such transformations. Since the angles and 
the sides (given by the moduli of the elements of the mixing matrix) in these triangles 
remain unchanged, they are phase convention independent and are physical observables. 
Consequently they can be measured directly in suitable experiments. The area of the unitarity 
triangle is related to the measure of CP violation Jcp [0, ll^ : 



I Jcp 1= 2 • Aa, 

where A^. denotes the area of the unitarity triangle. 

The construction of the unitarity triangle proceeds as follows: 



[1.26] 



We note first that 



-AX^ + 0{X'^). 

Thus to an excellent accuracy V^^V^^ is real with IV^^V^^^I = ^A^. 



Keeping 0{X^) corrections and rescaling all terms in ( |1.25 ) by AX"^ we find 



AX"^ 

with Q and fj defined in ( |1.21 ). 



AX^ 



VM = l-{g + ifj) 



(1.27) 



(1.28) 



Thus we can represent ( 1.25 ) as the unitarity triangle in the complex {Q,r]) plane as 
shown in fig. 



A=(p,r|) 



p^ri 




C=(0,0) 



l-p4r| 



B=(1,0) 



Figure 1: Unitarity Triangle. 



Let us collect useful formulae related to this triangle: 
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Using simple trigonometry one can express sin(2(/)j), 0j = a,/3,7, in terms of ??) as 
follows: 



sin(2a) 



sin(2/3) 
sin(27) = 



(^2 +r?2)((l - +f/2) 
2??(1 - Q) 



- (1 _ g)2 + 

2^7/ 2Qrj 



+ rf + r/2 



:i.29) 
:i.3o) 
:i.3i) 



The lengths CA and BA in the rescaled triangle to be denoted by ii^ and Rt, respec- 
tively, are given by 



Rb 



Vub 



Vcb 



l-^)2+r?2 



A 



Vtd 



Yr 



Cb 



[1.32] 



:i.33) 



The angles /3 and 7 of the unitarity triangle are related directly to the complex phases 
of the CKM-elements Vtd and Vub, respectively, through 



V 



td 



\ytd\e-'^ 



V 



ub 



\Vub\e-'^. 



:i.34) 



The angle a can be obtained through the relation 

a + /3 + 7 = 180° 



expressing the unitarity of the CKM-matrix. 



:i.35) 



The triangle depicted in fig. || together with \ Vus\ and \ Vcb\ gives a full description of the 
CKM matrix. Looking at the expressions for Rb and Rt, we observe that within the Standard 
Model the measurements of four CP conserving decays sensitive to | Vus |, | V^b |, | Vcb \ and 
I Vtd I can tell us whether CP violation (ry 7^ 0) is predicted in the Standard Model. This is a 
very remarkable property of the Kobayashi-Maskawa picture of CP violation: quark mixing 
and CP violation are closely related to each other. 



1.4 Grand Picture 

What do we know about the CKM matrix and the unitarity triangle on the basis of tree level 
decays? A detailed answer to this question can be found in the reports of the Particle Data 



Group 1^] as well as other reviews [14, 15], where references to the relevant experiments and 
related theoretical work can be found. In particular we have 



IK. 



A = 0.2205 ± 0.0018 



\Vcb\ = 0.040 ±0.002, 



(1.36) 
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\Vcb\ 

Using (|rT| and ( [1321) we find then 



- 0.089 ±0.016, iKbl = (3.56 ±0.56) • 10"^. (1.37) 



^ = 0.826 ±0.041, i?;, = 0.39 ± 0.07 . (1.38) 

This tells us only that the apex A of the unitarity triangle lies in the band shown in fig. 

In order to answer the question where the apex A lies on this "unitarity clock" we have 
to look at different decays. Most promising in this respect are the so-called "loop induced" 
decays and transitions which are the subject of several sections in these lectures and CP 
asymmetries in B-decays which will be briefly discussed in Section 8. These two different 



0.5 



Tj 



-0.5 



-0.5 0.5 

P 



Figure 2: "Unitarity Clock". 

routes for explorations of the CKM matrix and of the related unitarity triangle may answer 
the important question, whether the Kobayashi-Maskawa picture of CP violation is correct 
and more generally whether the Standard Model offers a correct description of weak decays 
of hadrons. Indeed, in order to answer these important questions it is essential to calculate 
as many branching ratios as possible, measure them experimentally and check if they all can 
be described by the same set of the parameters (A,^, £», ??). In the language of the unitarity 
triangle this means that the various curves in the (^, ff) plane extracted from different decays 
should cross each other at a single point as shown in fig. |3|. Moreover the angles (a, /3, 7) 
in the resulting triangle should agree with those extracted one day from CP-asymmetries in 
B-decays. For artistic reasons the value of fj in fig. ^ has been chosen to be higher than the 
fitted central value f) 0.35. 
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-1 ^0 1 2 

P 

Figure 3: The ideal Unitarity Triangle. 

On the other hand if new physics contributes to weak decays the different curves based 
on the Standard Model expressions, will not cross each other at a single point and the angles 
(a, P, 7) extracted one day from CP-asymmetries in B-decays will disagree with the ones 
determined from rare K and B decays. Clearly the plot in fig. ^ is highly idealized because 
in order to extract such nice curves from various decays one needs perfect experiments and 
perfect theory. One of the goals of these lectures is to identify those decays for which at least 
the theory is under control. For such decays, if they can be measured with a sufficient preci- 
sion, the curves in fig. ^ are not fully unrealistic. Let us then briefly discuss the theoretical 
framework for weak decays. 

2 Theoretical Framework 

2.1 OPE and Renormalization Group 

The basis for any serious phenomenology of weak decays of hadrons is the Operator Product 
Expansion (OPE) which allows to write the effective weak Hamiltonian simply as 

follows 

^e// = ^E^CKMa(/x)g. . (2.1) 

Here G_p is the Fermi constant and Qi are the relevant local operators which govern the 
decays in question. They are built out of quark and lepton fields. The Cabibbo-Kobayashi- 
Maskawa factors V^^i^jv/ @' HI ^^'^ Wilson coefficients Ci jl^ describe the strength with 
which a given operator enters the Hamiltonian. An amplitude for a decay of a given meson 
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M = K, B, .. into a final state F = ttuv^ tttt, DK is then simply given by 

A{M ^ F) = {F\neff\M) = ^Y.^^KMCi{lj){F\QmM), (2.2) 

where {F\Qi{^)\M) are the hadronic matrix elements of Qi between M and F. 

The essential virtue of OPE is this one. It allows to separate the problem of calculating the 
amplitude A{M F) into two distinct parts: the short distance (perturbative) calculation 
of the coefficients Ci{p,) and the long-distance (generally non-perturbative) calculation of 
the matrix elements {Qi{fi)). The scale fi separates the physics contributions into short 
distance contributions contained in Ci{fi) and the long distance contributions contained in 
{Qi{li)). Thus Ci include the top quark contributions and contributions from other heavy 
particles such as W, Z-bosons and charged Higgs particles or supersymmetric particles in the 
supersymmetric extensions of the Standard Model. Consequently Ci{^) depend generally on 
TTT-f and also on the masses of new particles if extensions of the Standard Model are considered. 
This dependence can be found by evaluating so-called hox and penguin diagrams with full 
W-, Z-, top- and new particles exchanges and properly including short distance QCD effects. 
The latter govern the //-dependence of Ci{fi). 

The value of fi can be chosen arbitrarily but the final result must be independent. 
Therefore the //-dependence of Cj(/i) has to cancel the //-dependence of (Qi (/*))• In other 
words it is a matter of choice what exactly belongs to Cj(//) and what to (Qi(//)). This 
cancellation of //-dependence involves generally several terms in the expansion in ( |2.2| ). The 
coefficients Ci(//) depend also on the renormalization scheme. This scheme dependence must 
also be cancelled by the one of {Qi{f^)) so that the physical amplitudes are renormalization 
scheme independent. Again, as in the case of the //-dependence, the cancellation of the 
renormalization scheme dependence involves generally several terms in the expansion ( [2.21 ). 

Although /i is in principle arbitrary, it is customary to choose /i to be of the order 
of the mass of the decaying hadron. This is 0{mh) and 0{mc) for B-decays and D-decays 
respectively. In the case of K-decays the typical choice is// = 0(1 — 2 GeV) instead oiO{mK), 
which is much too low for any perturbative calculation of the couplings Ci . Now due to the 
fact that // <C Mw,z, fnt, large logarithms lnMw//t compensate in the evaluation of Ci(/t) the 
smallness of the QCD coupling constant as and terms a"(lnMw//i)"', a"(lnMw//i)"'~^ etc. 
have to be resummed to all orders in before a reliable result for Cj can be obtained. This 
can be done very efficiently by means of the renormalization group methods. The resulting 
renormalization group improved perturbative expansion for Ci{jjL) in terms of the effective 
coupling constant ots{ii) does not involve large logarithms and is more reliable. 

All this looks rather formal but in fact should be familiar. Indeed, in the simplest case of 
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the /3-decay, Ti-eff takes the famihar form 



ni% = ^ cose,[^i7M(l - 75)d 67^(1 - 75)z^e] , (2.3) 



where Vud has been expressed in terms of the Cabibbo angle. In this particular case the Wilson 
coefficient is equal unity and the local operator, the object between the square brackets, is 



given by a product of two V — A currents. Equation (2.3) represents the Fermi theory for 
/3-decays as formulated by Sudarshan and Marshak [|l8| and Feynman and Gell-Mann |1S] 
forty years ago, except that in (|2.3|) the quark language has been used and following Cabibbo 
a small departure of Vud from unity has been incorporated. In this context the basic formula 



(2.1) can be regarded as a generalization of the Fermi Theory to include all known quarks and 
leptons as well as their strong and electroweak interactions as summarized by the Standard 
Model. 

Due to the interplay of electroweak and strong interactions the structure of the local 
operators is much richer than in the case of the /3-decay. They can be classified with respect 
to the Dirac structure, colour structure and the type of quarks and leptons relevant for a 
given decay. Of particular interest are the operators involving quarks only. In the case of the 
AS = 1 transitions the relevant set of operators is given as follows: 

Current— Current : 

Ql = {SaU/3)v~A {upda)v-A Q2 = {su)v-A {ud)v-A (2.4) 

QCD Penguins : 

Q3 = {sd)v-A X! i^^)v-A Qi = {sadp)v-A ^ {qpqa)v-A (2.5) 

q=u,d,s q=u,d,s 

Qb = {sd)v-A i9q)v+A Qe = {Sadf^)v-A i^p1a)v+A (2.6) 

q=u,d,s q=u,d,s 

Electroweak Penguins : 

3 3 

= 2 {^d)v-A ^9 iQl)v+A 2 (iq{Wa)v+A (2.7) 

q=u,d,s q=u,d,s 

3 3 

Q9 = ^{sd)v-A Y eg(gg)y-A QlO = ^ {Sadf3)v-A Y {q/sqaW-A ■ (2.8) 
q=u,d,s q=u,d,s 

Here, Cg denotes the electric quark charges reflecting the electroweak origin of Qy, . . . , Qiq. 

Clearly, in order to calculate the amplitude A(M F), the matrix elements {Qi{^)) 
have to be evaluated. Since they involve long distance contributions one is forced in this case 
to use non-perturbative methods such as lattice calculations, the 1/N expansion (N is the 
number of colours), QCD sum rules, hadronic sum rules, chiral perturbation theory and so 
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on. In the case of certain B-meson decays, the Heavy Quark Effective Theory (HQET) also 
turns out to be a useful tool. Needless to say, all these non-perturbative methods have some 
limitations. Consequently the dominant theoretical uncertainties in the decay amplitudes 
reside in the matrix elements {Qi{fi)). 

The fact that in most cases the matrix elements (Qi(/x)) cannot be reliably calculated at 
present, is very unfortunate. One of the main goals of the experimental studies of weak decays 
is the determination of the CKM factors Vckm and the search for the physics beyond the 
Standard Model. Without a reliable estimate of (QiifJ-)) this goal cannot be achieved unless 
these matrix elements can be determined experimentally or removed from the final measurable 
quantities by taking the ratios or suitable combinations of amplitudes or branching ratios. 
However, this can be achieved only in a handful of decays and generally one has to face 
directly the calculation of {Qi{ii)). We will discuss these issues later on. 



2.2 Inclusive Decays 

So far I have discussed only exclusive decays. It turns out that in the case of inclusive decays 
of heavy mesons, like B-mesons, things turn out to be easier. In an inclusive decay one sums 
over all (or over a special class) of accessible final states so that the amplitude for an inclusive 
decay takes the form: 

A{B^X) = ^Y1 V^^Ml^){fmi^)\B) . (2.9) 

At first sight things look as complicated as in the case of exclusive decays. It turns out, 
however, that the resulting branching ratio can be calculated in the expansion in inverse 
powers of mb with the leading term described by the spectator model in which the B-meson 
decay is modelled by the decay of the 6-quark: 

BiiB ^ X) =Br{b^ q) + Oi-^) . (2.10) 



This formula is known under the name of the Heavy Quark Expansion (HQE) ||2^. Since 
the leading term in this expansion represents the decay of the quark, it can be calculated in 
perturbation theory or more correctly in the renormalization group improved perturbation 
theory. It should be realized that also here the basic starting point is the effective Hamiltonian 



(2.1) and that the knowledge of Ci{^) is essential for the evaluation of the leading term in 
( 2.10| ). But there is an important difference relative to the exclusive case: the matrix elements 
of the operators Qi can be "effectively" evaluated in perturbation theory. This means, in 
particular, that their /x and renormalization scheme dependences can be evaluated and the 
cancellation of these dependences by those present in Cj(/i) can be investigated. 
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Clearly in order to complete the evaluation of Br{B X) also the remaining terms in 
( 2.10| ) have to be considered. These terms are of a non-perturbative origin, but fortunately 



they are suppressed by at least two powers of nib. They have been studied by several authors 
in the literature with the result that they affect various branching ratios by less then 10% 
and often by only a few percent. Consequently the inclusive decays give generally more 
precise theoretical predictions at present than the exclusive decays. On the other hand their 
measurements are harder. There are of course some important theoretical issues related to 
the validity of HQE in plc| ) which appear in the literature under the name of quark-hadron 



duality. Since these matters are rather involved I will not discuss them here. 

2.3 Status of NLO Calculations 

In order to achieve sufficient precision for the theoretical predictions it is desirable to have 
accurate values of Ci{ji). Indeed it has been realized at the end of the eighties that the 
leading term (LO) in the renormalization group improved perturbation theory, in which the 
terms Q"(lnMw/^)"' are summed, is generally insufficient and the inclusion of next-to-leading 
corrections (NLO) which correspond to summing the terms a" (luMw//^)""^ is necessary. In 
particular, unphysical left-over /i-dependences in the decay amplitudes and branching ratios 
resulting from the truncation of the perturbative series are considerably reduced by including 
NLO corrections. These corrections are known by now for the most important and interesting 
decays and will be taken into account in these lectures. The review of all existing NLO 
calculations can be found in [^, 

2.4 Penguin— Box Expansion 

The rare and CP violating decays of K and B mesons are governed by various penguin and box 
diagrams with internal top quark and charm quark exchanges. Some examples are shown in 
fig. Evaluating these diagrams one finds a set of basic universal (process independent) rrit- 
dependent functions Fr{xt) |22] where xt = /M^. ExpUcit expressions for these functions 
will be given below. 



It is useful to express the OPE formula ( p.2| ) directly in terms of the functions Fj.{xt) |2J 
A{M -^F) = Po{M -^F) + J2 Pr{M ^ F) Fr{xt), (2.11) 

r 

where the sum runs over all possible functions contributing to a given amplitude. Pq summa- 
rizes contributions stemming from internal quarks other than the top, in particular the charm 
quark. The coefficients Pq and Pr are process dependent and include QCD corrections. They 
depend also on hadronic matrix elements of local operators and the relevant CKM factors. I 
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u,c,t 



u,c,t 




Figure 4: Typical Penguin and Box Diagrams. 



would like to call ( 2.11| ) Penguin-Box Expansion (PBE). We will encounter many examples 
of PBE in the course of these lectures. 

Originally PBE was designed to expose the mt-dependence of FCNC processes [^3|. After 
the top quark mass has been measured precisely this role of PBE is less important. On the 
other hand, PBE is very well suited for the study of the extentions of the Standard Model 
in which new particles are exchanged in the loops. If there are no new local operators the 
mere change is to modify the functions Fr{xt) which now acquire the dependence on the 
masses of new particles such as charged Higgs particles and supersymmetric particles. The 
process dependent coefficients Pq and Pr remain unchanged. The effects of new physics can 
be then transparently seen. However, if new effective operators with different Dirac and 
colour structures are present the values of Pq and Pr are modified. 

Let us denote by Bq, Cq and Dq the functions Fr{xt) resulting from AF = 1 {F stands 
for fiavour) box diagram, -penguin and 7-penguin diagram respectively. These diagrams 



14 



are gauge dependent and it is useful to introduce gauge independent combinations ||2^ 

Xo = C7o-4So, Yo = Co-Bo, z^ = Co + \dq (2.12) 

Then the set of gauge independent basic functions which govern the FCNC processes in the 
Standard Model is given to a very good approximation as follows: 

S„(x,)=2A6 ( i^Q^y ) ■ S„{x,) = x, (2,13) 



The first three functions correspond to IS.F = 2 box diagrams with (t, i), (c, c) and (t, c) 
exchanges. -Eg results from QCD penguin diagram with off-shell gluon, D'q and Eq from 
7 and QCD penguins with on-shell photons and gluons respectively. The subscript "0" 
indicates that these functions do not include QCD corrections to the relevant penguin and 
box diagrams. 

In the range 150 GeV < m-t < 200 GeV these approximations reproduce the exact expres- 
sions to an accuracy better than 1%. These formulae will allow us to exhibit elegantly the 
TTit dependence of various branching ratios in the phenomenological sections of these lectures. 
Exact expressions for all functions can be found in Q. 

Generally, several basic functions contribute to a given decay, although decays exist which 
depend only on a single function. We have the following correspondence between the most 
interesting FCNC processes and the basic functions: 





— ^'^-mixin; 


y 


Soixt), So{Xc,Xt), So{Xc) 




— iJ'^-mixing 




Soixt) 


K - 


TTiyv', B 




Xo{xt) 


Kl 


fifi, B ^ 
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Yoixt) 


Kl 


7r^e~^e~ 




Yo{xt), Zo{xt), Eo{xt) 


e' 






Xo{xt), Yo{xt), Zo{xt), Eo{xt) 


B - 






D',{xt), E'^{xt) 


B - 






Yoixt), 2'o(xj), -Eo(xt), DQ{xt), EQ{xt 
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3 Particle- Antiparticle Mixing and CP Violation 



3.1 Preliminaries 



Let us next discuss particle-antiparticle mixing which in the past has been of fundamental 
importance in testing the Standard Model and often has proven to be an undefeatable chal- 
lenge for suggested extensions of this model. Let us just recall that from the calculation of 
the -RTl — Ks mass difference, Gaillard and Lee [^] were able to estimate the value of the 
charm quark mass before charm discovery. On the other hand — mixing ||2^ gave the 
first indication of a large top quark mass. Finally, particle-antiparticle mixing in the — K^ 
system offers within the Standard Model a plausible description of CP violation in Kj^ — > tttt 



discovered in 1964 [26|. 

In this section we will predominantly discuss the parameter e describing the indirect 
CP violation in the K system and the mass differences AM^ ,, which describe the size of 
B'^g — B^g mixings. In the Standard Model these phenomena appear first at the one-loop 
level and as such they are sensitive measures of the top quark couplings Vti{i = d,s,b) and 
and in particular of the phase 6 = j. They allow then to construct the unitarity triangle. 

Let us next enter some details. The following subsection borrows a lot from |27, 28]. A 



nice review of CP violation can also be found in |29]. 



s W d 
-< •^/\/\/\A/\/V\y# <— 



u,c,t 



u,c,t 



— *wvww\> — >- 
d W s 

(a) 



u,c,t 

— « — 



u,c,t 
(b) 



d 



W 



Figure 5: Box diagrams contributing to — mixing in the Standard Model. 



3.2 Express Review of — Mixing 

= (sd) and = (sd) are flavour eigenstates which in the Standard Model may mix via 
weak interactions through the box diagrams in fig. ||. We will choose the phase conventions 
so that 

CP\K^) = -\K^), CP\K^) = -\K^). (3.1) 
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In the absence of mixing the time evolution of \K^{t)) is given by 

= |if°(0)) exp(-im) , H = M-i^ , (3.2) 

where M is the mass and F the width of K^. Similar formula for exists. 

On the other hand, in the presence of flavour mixing the time evolution of the — 
system is described by 



where 

^ "Zll-Z^ 1V1X2 

2 \ Moi-i^ Moo-i^ 



21 - -(^22 - 

with M and f being hermitian matrices having positive (real) eigenvalues in analogy with M 
and r. Mij and are the transition matrix elements from virtual and physical intermediate 
states respectively. Using 

M21 = Mi*2 , F21 = , (hermiticity) (3.5) 

Mii=M22 = M, Fii=F22 = F, (CPT) (3.6) 



we have 



We can next diagonalize the system to find: 
Eigenstates: 



where £ is a small complex parameter given by 

1-e I Mi*2 - 4^2 AM-^AF ^ ^ ^ ^ 
= ^2 ^_J2 ^ 2^ = rexp(iAv) . (3.9) 

1 + e ^Mi2-iiFi2 2Mi2-zFi2 ^ ^ ^ 

with AF and AM given below. 
Eigenvalues: 

ML,5 = M±Reg Fi,5 = rT2ImQ (3.10) 
where 



Q = ^(Mi2 - 4ri2)(Mf2 - i^Tl,). (3.11) 
Consequently we have 

AM = Ml - Ms = 2ReQ AF = F^ - F5 = -4ImQ. (3.12) 
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It should be noted that the mass eigenstates Ks and differ from CP eigenstates 

Ki = - CP\Ki) = \K{) , (3.13) 

K2 = -^(if° + K°), CP\K2) = -\K2) , (3.14) 

by a small admixture of the other CP eigenstate: 



Ki + eK2 K2 + eKi 



It should be stressed that the small parameter e depends on the phase convention chosen 
for and K^. Therefore it may not be taken as a physical measure of CP violation. On 
the other hand Ree and r are independent of phase conventions. In particular the departure 
of r from 1 measures CP violation in the K° - K° mixing: 

r = l+ „ Im ( ^) . (3.16) 



4|Mi2|2 + |ri2|2 Vri2 

Since e is 0(10^^), one has to a very good approximation: 

AMk = 2ReMi2, ATk = 2Reri2 , (3.17) 

where we have introduced the subscript K to stress that these formulae apply only to the 
K° - K° system. 

The -f^L — Ks mass difference is experimentally measured to be |8| 

AMk = M{Ki^) - M{Ks) = (3.489 ± 0.009) • 10^^^ GeV . (3.18) 

In the Standard Model roughly 70% of the measured AMk is described by the real parts of 
the box diagrams with charm quark and top quark exchanges, whereby the contribution of 
the charm exchanges is by far dominant. This is related to the smallness of the real parts 
of the CKM top quark couplings compared with the corresponding charm quark couplings. 
Some non-negligible contribution comes from the box diagrams with simultaneous charm 
and top exchanges. The remaining 20% of the measured AMk is attributed to long distance 



contributions which are difficult to estimate [30|. Further information with the relevant 
references can be found in |31]. 

The situation with AP^^ is rather different. It is fully dominated by long distance effects. 
Experimentally one has AP^ ~ — 2AMi^. 
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3.3 The First Look at e and s' 



Since a two pion final state is CP even while a three pion final state is CP odd, and Ki, 
preferably decay to 27r and Svr, respectively via the following CP conserving decay modes: 



i^L 37r (via Ka), Ks ^ 27r (via Ki). 



(3.19) 



This difference is responsible for the large disparity in their life-times. A factor of 579. 
However, Ki^ and are not CP eigenstates and may decay with small branching fractions 
as follows: 

ETl 27r (via Ki), Ks Svr (via Ka). (3.20) 

This violation of CP is called indirect as it proceeds not via explicit breaking of the CP 
symmetry in the decay itself but via the admixture of the CP state with opposite CP parity 
to the dominant one. The measure for this indirect CP violation is defined as 

A{Ki, (7r7r)/=o) 



A{Ks (vr^)/=o) ' 
Following the derivation in [^] one finds 

exp(i7r/4) , 

' ' (ImMi2 + 2eReMi2) , 



Imyln 



(3.21) 



(3.22) 



^/2AMi^ 

where the term involving ImMia represents e defined in (| 
dence of the term involving ^ cancells the convention dependence of e so that e is free from 
this dependence. 



Revlo 

The phase convention depen- 



Ki cX K 



z K 



1 



indirect: z 



direct : e' 



TCTI 



TTTC 



Figure 6: Indirect versus direct CP violation in Kl — > vrvr. 

While indirect CP violation reflects the fact that the mass eigenstates are not CP eigen- 
states, so-called direct CP violation is realized via a direct transition of a CP odd to a CP 
even state or vice versa (see fig. ^). A measure of such a direct CP violation in Kl — > tttt is 
characterized by a complex parameter e' defined as 



i=In, (^J exp(z$,0, 



vr 



+ ^2 



^0, 



(3.23) 
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where the isospin amphtudes Aj in K ^ vrvr decays are introduced through 

A{K+ ^ vr+vrO) = ^1^26^^^ (3.24) 

A{K^ ^ vr+vr-) = /|^oe''° + V^^2e^'^ (3.25) 

^(^0 ^ vrOvrO) = ^l^oe'^" - 2^A2e'^' . (3.26) 

Here the subscript / = 0,2 denotes states with isospin 0,2 equivalent to A/ = 1/2 and 
A/ = 3/2 transitions, respectively, and 5q^2 are the corresponding strong phases. The weak 
CKM phases are contained in ^'iid A2. The isospin amplitudes Aj are complex quantities 
which depend on phase conventions. On the other hand, e' measures the difference between 
the phases of A2 and Aq and is a physical quantity. The strong phases ^0,2 can be extracted 
from TTvr scattering. Then ~ 7r/4. 

Experimentally e and e' can be found by measuring the ratios 

_ A{Ki^ ^ ^"vrO) _ A{Ki^ ^ ^+^-) 

Indeed, assuming e and e' to be small numbers one finds 

2e' e' 

r]oo = e- ^ ~ e - 2e', ry+_ = e + — = ~ e + e' (3.28) 

i — yu; 1 + uj/\/ 2 

where experimentally uj = ReA2/ReAo = 0.045. 

In the absence of direct CP violation 7700 = The ratio e' /e can then be measured 

through 



Voo 



V+- 



-1-6 Re(^) . (3.29) 



3.4 Basic Formula for e 



With all this information at hand let us derive a formula for e which can be efficiently used in 
pheneomenological applications. The off-diagonal element M12 in the neutral i^-meson mass 
matrix representing K^-K^ mixing is given by 

2mKMl2 = {K^lHesi^S = 2)\K^) , (3.30) 

where Wgff (A^ = 2) is the effective Hamiltonian for the AS" = 2 transitions. That M*2 and 
not M12 stands on the l.h.s of this formula, is evident from (3.7). The factor 2mx reflects 
our normalization of external states. 
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To lowest order in electroweak interactions AS" = 2 transitions are induced through the 
box diagrams of fig. Including QCD corrections one has 



n_/AS=2 



167r2 ^ 



X^r]iSo{xc) + Xi.r]2So{xt) + 2XcXtr]3So{xc, xt] 



-2/9 



J- H : -'3 



in 



Q{AS = 2) + h.c. 



(3.31) 



where Aj = V*gV^^, fi < fic = 0{mc) and is the strong coupling constant in an eff'ective 



three flavour theory. In (3.31), the relevant operator 

Q(A5 = 2) = {sd)v-A{sd)v-A, 



(3.32) 



is multiplied by the corresponding coefficient function. This function is decomposed into a 
charm-, a top- and a mixed charm-top contribution. The functions are given in (|2l^)and 
(EH)- 

Short-distance QCD effects are described through the correction factors r/i, r/2, % and the 
explicitly a^-dependent terms in ( 3.31| ). The NLO values of rn are given as follows [31, 33]: 



r?i = 1.38 ±0.20, 772 = 0.57 ±0.01, 



r?3 = 0.47 ± 0.04 . 



(3.33) 



The quoted errors reflect the remaining theoretical uncertainties due to leftover /^-dependences 
at 0{a^g) and ^jjg, the scale in the QCD running coupling. 

Deflning the renormalization group invariant parameter Bk by 



Bk = BM af\p) 



-2/9 



{K\sd)v-A{sd)v^A\K'' 
and using ( 3.31| ) one flnds 

M12 



-BMFlml 



F^BkuikM^ \xfmSo{xc) + Xfv2So{xt) + 2A*A*%So(x„ x*) 



127r2^ 



(3.34) 
(3.35) 

(3.36) 



where is the i^-meson decay constant and rriK the ET-meson mass. 

To proceed further we neglect the last term in ( ^.22 ) as it constitutes at most a 2% 
correction to e. This is justified in view of other uncertainties, in particular those connected 
with Bk- Inserting (|3^ ) into (|]2|) we find 



e = CeBK^mXt{ReXc [r]iSo{xc) - mSoixcXt)] - ReXtmSoixt)} exp{i7r / 4) 



(3.37) 



where we have used the unitarity relation ImA* = IniA^ and have neglected ReAf/ReAc 
O(A^) in evaluating Im(A*AJ'). The numerical constant Cg is given by 



C. = ^i#=^ =3.84. 10^ 



(3.38) 
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To this end we have used the experimental value of AMj^ in ( 3.1q ). 

Using the standard parametrization of ( |1.5| ) to evaluate ImAj and ReA,, setting the values 
for si2, si3, S23 and mt in accordance with experiment and taking a value for Bk (see below), 
one can determine the phase 6 by comparing ( |3.37| ) with the experimental value for e 

Eexp = (2.280 ± 0.013) • 10"^ expi^s, = J. (3.39) 

Once 6 has been determined in this manner one can find the apex (^,77) of the unitarity 
triangle in fig. |l] by using 

cos 5, r] = — sin 6 (3.40) 



and 



S12S23 S12S23 



^^(l-Y), r/ = ry(l-^). (3.41) 



For a given set {su, S13, S23, mt, Bk) there are two solutions for 5 and consequently two 
solutions for {g,fi)- This will be evident from the analysis of the unitarity triangle discussed 
in detail below. 

Finally we have to say a few words about the non-perturbative parameter Bj^, the main 
uncertainty in this analysis. Reviews are given in |3^, Here we only collect in table |^ 
values for Bk obtained in various non-perturbative approaches. In our numerical analysis 
presented below we will use 

= 0.80 ±0.15 (3.42) 
which is in the ball park of various lattice and large-N estimates. 

3.5 Basic Formula for B^-B^ Mixing 

The strength of the B'^ ^ — B^ ^ mixings is described by the mass differences 

AMd,s = M^' - Mf^'' (3.43) 

with "H" and "L" denoting Heavy and Light respectively. In contrast to AMk, in this 
case the long distance contributions are estimated to be very small and AAI^^s is very well 
approximated by the relevant box diagrams. Moreover, due rriu^c ^ mt only the top sector 
can contribute significantly to B^^ — B^ ^ mixings. The charm sector and the mixed top-charm 
contributions are entirely negligible. 

AMfi^g can be expressed in terms of the off-diagonal element in the neutral B-meson mass 
matrix by using the formulae developed previously for the K-meson system. One finds 

AMq = 2\M[^\, q = d,s. (3.44) 
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Table 1: obtained using various methods. WA stands for recent world avarage. 



Method 


Br 


Reference 


Chiral QM 


1.1 ±0.2 


m 


Lattice (APE) 


0.93 ±0.16 


[|7 


] 


Lattice (JLQCD) 


0.86 ± 0.06 




Lattice (GKS) 


0.85 ± 0.05 




Lattice (WA) 


0.89 ±0.13 




i] 


Large-N 


0.70 ±0.10 


|41 




1 


Large-N 


0.4-0.7 


1 


4|] 


QCDS 


0.5 - 0.6 


1 


44 


] 


CHPTH 


0.42 ± 0.06 


1 


m 


QCD HD 


0.39 ±0.10 


[ 




SU(3)+PCAC 


0.33 


1 


47 


] 



This formula differs from AMk = 2ReMi2 because in the B-system <^ Mi2- 

The off-diagonal term M12 in the neutral S-meson mass matrix is then given by a formula 
analogous to (IH) 

(3.45) 



2mBjMi$\ = \{B^^\n,s{AB = 2)\B% 



where in the case of B^ — B^ mixing 



n,AB = 2 



167r2" 



-6/23 



Q{AB = 2) + h.c. 



Here ^b = ^{n^b), 



QiAB = 2) = {bd)v-A{bd)v-A 



and |3|1 



(3.46) 



(3.47) 



(3.48) 



r]B = 0.55 ±0.01. 

In the case of B^ — B^ mixing one should simply replace d ^ s \n (3.46) and ( p. 47 ) with all 
other quantities unchanged. 

Defining the renormalization group invariant parameters I3q in analogy to ( 3.34| ) and 
( 3.35| ) one finds using (3.46) 



AM, = ^^j,mB,{BB,Fl^)M^So{xt)\Vtq\\ 



(3.49) 
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where Fb^ is the i?q-meson decay constant. This imphes two useful formulae 



AMd = 0.50/ps • 



200 MeV 



"T-t(^t) 

170 GeV 



1.52 



td\ 



8.8 • 10- 



0.55 



(3.50) 



and 



AM^ = 15.1/ps 



BbFb, 



240 MeV 



mt[mt) 



170 GeV 



1.52 r 



Wt 



n 2 r 



ts\ 



0.040 



0.55 



(3.51) 



There is a vast literature on the calculations of Fb^ and Bd. The most recent lattice 
results are summarized in [^^. They are compatible with the results obtained with the help 
of QCD sum rules 1^^. In our numerical analysis we will use the value for Fb^^/Bb^ given 
in table 0. The experimental situation on AMg^ g is also given there. 



4 Standard Analysis of the Unitarity Triangle 

4.1 Basic Procedure 

With all these formulae at hand we can now summarize the standard analysis of the unitarity 
triangle in fig. |l[ It proceeds in five steps. 
Step 1: 

From 6 — > c transition in inclusive and exclusive leading B meson decays one finds \Vcb\ 
and consequently the scale of the unitarity triangle: 



cb\ 



MVcb\ = X^A 



(4.52) 



Step 2: 

From b ^ u transition in inclusive and exclusive B meson decays one finds |l^fe/V^b| and 
consequently the side CA = Rf, of the unitarity triangle: 



Vub 



Vcb 



Rf^ = J -\- fj2 = 4.44 . 



Vub 



Vcb 



(4.53) 



Step 3: 



From the experimental value of e ( 3.39| ) and the formula ( p. 37 ) one derives, using the 
approximations (|1.22| )-( 1.24 ), the constraint 



(1 - g)A^rj2Soixt) + Po(e) A^Bk = 0.224, 



(4.54) 



where 



(4.55) 



w 
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Po{£) = 0.31 lb 0.05 summarizes the contributions of box diagrams with two charm quark 
exchanges and the mixed charm-top exchanges. The main uncertainties in the constraint 
(4.54) reside in Bk and to some extent in yl^ which multiphes the leading term. Equation 
( 4.54 ) specifies a hyperbola in the (^, ?]) plane. This hyperbola intersects the circle found in 
step 2 in two points which correspond to the two solutions for 6 mentioned earlier. This is 
illustrated in fig. ^. The position of the hyperbola (4.54) in the (^, r/) plane depends on rrit, 
\Vcb\ = AX^ and Bk- With decreasing mt, \Vcb\ and IBk the e- hyperbola moves away from 
the origin of the {g, fj) plane. 




Figure 7: Schematic determination of the Unitarity Triangle. 



Step 4: From the observed B^ — B^ mixing parametrized by AMd the side BA = Rt of 
the unitarity triangle can be determined: 



1 W, 



td\ 



A Wr 



1.0 



cb\ 



td\ 



i.8 ■ 10-3 



0.040 



IK 



with 



td\ 



.8 • 10" 



200 MeV 



170 GeV 



0.76 



cb\ 



AM, 



0.50/ps 



0.5 



/0.55 
r]B 



(4.56) 



(4.57) 



mt(mt) 

Since rrii, AM^ and rjB are already rather precisely known, the main uncertainty in the 
determination of \ Vtd\ from B^ — mixing comes from Fb^^J Bb^- Note that Rt suffers from 
additional uncertainty in \Vcb\, which is absent in the determination of \Vtd\ this way. The 
constraint in the {g, ff) plane coming from this step is illustrated in fig. |7[ 

Step 5: 
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Table 2: Collection of input parameters. 



Quantity 



Central 



Error 



Reference 



\Vuh\ 

Bk 
mt 
AMs 



0.040 
3.56 • 10-3 

0.80 
200 MeV 
165 GeV 
0.471 ps-i 
> 12.4 ps-i 
1.14 



±0.002 
±0.56 • 10-3 
±0.15 
±40 MeV 
±5 GeV 
±0.016 ps-i 
95%C.L. 
±0.08 



in 

See Text 

m 
m 
m 
m 

mm 



The measurement of — mixing parametrized by AMg together with AM^ allows to 
determine Rt in a different way. Using ( p.49| ) one finds 



Using next AM™^^ = 0.487/ps and | Vts/Kbl™'^^ = 0.991 one finds a useful approximate 
formula 

One should note that mt and \Vcb\ dependences have been eliminated this way and that 
^ should in principle contain much smaller theoretical uncertainties than the hadronic ma- 
trix elements in AM^ and AM, separately. The most recent values relevant for ( [4.59| ) are 
summarized in table 0. 



4.2 Numerical Results 

4.2.1 Input Parameters 

The input parameters needed to perform the standard analysis of the unitarity triangle are 
given in table where mt refers to the running current top quark mass defined at /i = mf 
It corresponds to mf"'"" = 174.3 ± 5.1 GeV measured by CDF and DO @. 

4.2.2 Output of the Standard Analysis 



In what follows we will present two types of numerical analyses [35, p3[: 
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Figure 8: Unitarity Triangle 1999. 

• Method 1: The experimentally measured numbers are used with Gaussian errors and 
for the theoretical input parameters we take a flat distribution in the ranges given in 
table 

• Method 2: Both the experimentally measured numbers and the theoretical input pa- 
rameters are scanned independently within the ranges given in table ^. 

The results are shown in table |3|. The allowed region for f\) is presented in fig. ^ It is 
the shaded area on the right hand side of the solid circle which represents the upper bound 
for (AM)^/(AM)s. The hyperbolas give the constraint from e and the two circles centered 
at (0,0) the constraint from l^^t/lcbl- The white areas between the lower e-hyperbola and 
the shaded region are excluded by i?^ — mixing. We observe that the region ^ < is 
practically excluded. The results in fig. |^ correspond to a simple independent scanning of 
all parameters within one standard deviation. We find that whereas the angle (5 is rather 
constrained, the uncertainties in a and 7 are substantially larger: 



66° < a < 113° , 17° < /3 < 29° , 44° < 7 < 97° 



(4.60) 



The result for sin 2/3 is consistent with the recent measurement of CP asymmetry in B — > "0^5 



by CDF 1 54], although the large experimental error precludes any definite conclusion. 



Other studies of the unitarity triangle can be found in [14, 15, 55 
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Table 3: Output of the Standard Analysis. At = V*gVtd- 



Quantity 


Scanning 


Gaussian 


1 Vtd 1 /10-3 


7.0-9.8 


8.1 ±0.6 


1 Vts/Vd, 1 


0.975 - 0.991 


0.984 ± 0.004 


1 Vtd/Vts 1 


0.17-0.24 


0.201 ±0.017 


sin(2/3) 


0.57 - 0.84 


0.73 ±0.09 


sin(2Q;) 


-0.72 - 1.0 


-0.03 ±0.36 


sin(7) 


0.69 - 1.0 


0.89 ± 0.08 


ImAi/lO"^ 


1.04 - 1.63 


1.33 ±0.14 



4.3 Final Remarks 

In this section we have completed the determination of the CKM matrix. It is given by the 
values of \Vus\, \Vcb\ and \Vub\ in ( |1.36| ) and ( 1.37] ), the results in table |3| and the unitarity 



triangle shown in fig. Clearly the accuracy of this determination is not yet impressive. We 
should stress, however, that in a few years from now the standard analysis may give much 
more accurate results. In particular a single precise measurement of AMg will have a very 
important impact on the allowed area in the (g, f)) plane. Such a measurement should come 
from SLD and later from LHC. 

Having the values of CKM parameters at hand, we can use them to predict various 
branching ratios for rare and CP-violating decays. This we will do in the subsequent sections. 

5 e' /e in the Standard Model 
5.1 Preliminaries 

Direct CP violation remains one of the important targets of contemporary particle physics. 
In the case oi K ^ vrvr, a non-vanishing value of the ratio Re(e'/e) defined in ( p. 23 ) would 



give the first signal for direct CP violation ruling out superweak models [p7|]. Until recently 
the experimental situation on e'/e was rather unclear: 

R.(.7.) = ^ P3±7).10- (NA31)[581 
^ (7.4 ±5.9) • 10-4 (E731) ||. 



While the result of the NA31 collaboration at CERN [^] clearly indicated direct CP violation, 
the value of E731 at Fermilab [B^, was compatible with superweak theories [^] in which 
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e' je = 0. This controversy is now settled with the very recent measurement by KTeV at 
Fermilab [^] 

Re(e7e) = (28.0 ± 4.1) • 10"^ (KTeV) (5.2) 

which together with the NA31 result confidently establishes direct CP violation in nature. 
The grand average including NA31, E731 and KTeV results reads [^0| 

Re(e7e) = (21.8 ± 3.0) • 10""^ (5.3) 

very close to the NA31 result but with a smaller error. The error should be further reduced 
once the first data from NA48 collaboration at CERN are available and complete data from 
both collaborations have been analyzed. It is also of great interest to see what value for e' /e 
will be measured by KLOE at Frascati, which uses a different experimental technique than 
KTeV and NA48. 

There is a long history of calculations of e'/e in the Standard Model. The first calculation 
of e'/e for rn-t ^ Mw without the inclusion of renormalization group effects can be found in 
|l6l|| . Renormalization group effects in the leading logarithmic approximation have been first 
presented in [B^. For mt <C Mw only QCD penguins play a substantial role. First extensive 



phenomenological analyses in this approximation can be found in |63|. Over the eighties 
these calculations were refined through the inclusion of QED penguin effects for mt <C Mw 
p3, Bsl, 0, the inclusion of isospin breaking in the quark masses p5|, and through 



improved estimates of hadronic matrix elements in the framework of the 1/A^ approach |65]. 
This era of e'/e culminated in the analyses in ||6^, 70|, where QCD penguins, electroweak 
penguins (7 and penguins) and the relevant box diagrams were included for arbitrary top 
quark masses. The strong cancellation between QCD penguins and electroweak penguins for 
mt > 150 GeV found in these papers was confirmed by other authors [^]. 

During the nineties considerable progress has been made by calculating complete NLO 
corrections to e' |72|-|7^]. Together with the NLO corrections to e and — mixing |31 



32, 33|, this allowed a complete NLO analysis of e'/e including constraints from the observed 
indirect CP violation (e) and B^ ^ — B^ ^ mixings (AMd^s)- The improved determination of the 

and Vcb elements of the CKM matrix, the improved estimates of hadronic matrix elements 
using the lattice approach as well as other non-perturbative approaches and in particular the 
determination of the top quark mass mt had of course also an important impact on e'/e. 

Now, e'/e is given by ( p.5\ ) where in a crude approximation (not to be used for any serious 
analysis) 



Fe' « 13 



110 MeV "'2 



ms(2 GeV) 



R^^/^Vl O ,^0 4 0(3/2) / rnt ' ^'^ 
B, (l-a,+,0-0.4-i?8 (y65G^ 



MS 



340 MeV 



(5.4) 
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Here ^Irj+ri' ~ 0.25 represents isospin breaking corrections and Bi are hadronic parameters 
which we will define later on. This formula exhibits very clearly the dominant uncertainties 
in F;,! which reside in the values of nig, Bq'^'^\ B^^'^\ ^77+17' • Moreover, the partial 



cancellation between QCD penguin {Bq'^'^^) and electroweak penguin {B^^'^') contributions 



(3/2), 



requires accurate values of B^^"^^ and B^^'^' for an acceptable estimate of e'/e. Because of 
the accurate value mt(mt) = 165 ± 5 GeV, the uncertainty in e'/e due to the top quark mass 
amounts only to a few percent. A more accurate formula for Ff,i will be given below. 

Now, it has been known for some time that for central values of the input parameters the 
size of e'/e in the Standard Model is well below the NA31 value of (23.0 it 6.5) • 10"'^. Indeed, 



,(3/2) 



extensive NLO analyses with lattice and large~N estimates of B^ 
performed first in 



(1/2) 



1 and B, 



(3/2) 



1 



and after the top discovery in |77, 78, 79] have found e'/e in the 



ball park of (3 — 7) • 10 ^ for ms(2 GeV) ~ 130 MeV. On the other hand it has been stressed 
repeatedly in ||^, |7^ that for extreme values of 5, 



(1/2) „(3/2) 
1 Oo 



and rris still consistent with 



lattice, QCD sum rules and large-N estimates as well as sufficiently high values of ImAf and 
a|^, a ratio e'/e as high as (2 — 3) • 10~^ could be obtained within the Standard Model. Yet, it 
has also been admitted that such simultaneously extreme values of all input parameters and 
consequently values of e'/e close to the NA31 result are rather improbable in the Standard 
Model. Different conclusions have been reached in |8C], where values (1 — 2) • 10"'^ for e'/e can 
be found. Also the Trieste group 



, which calculated the parameters B^^'^'^ and B^^'^^ in 



the chiral quark model, found e'/e = (1.7 it 1.4) • 10^'^. On the other hand using an effective 
chiral lagrangian approach, the authors in |8^] found e'/e consistent with zero. 

After these general remarks let us discuss e'/e in explicit terms. Other reviews of e'/e 



can be found in |82, pl|. 



5.2 Basic Formulae 



The parameter e' is given in terms of the isospin amplitudes Aj in ( 3.23| ). Applying OPE to 
these amplitudes one finds 



where 



with 



Fe' ^ 

p(l/2) 
p(3/2) 



ImXt ■ F^,, 



p(l/2) _ p(3/2) 



exp(z$). 



(5.5) 
(5.6) 

(5.7) 
(5.8) 
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Here 



Since 



2|e|Reylo 



{{T:Ti)imK) 



Re^2 
Re^n' 



0, 



(5.9) 
(5.10) 



Ff,i and e' /e are real to an excellent approximation. The operators Qi have been given already 
in (p.4D-( p^ . The Wilson coefficient functions 2/i(/i) were calculated including the complete 
next-to-leading order (NLO) corrections in [|72|-||7^. The details of these calculations can 
be found there and in the review |||]. Their numerical values for corresponding to 



table 



0.119 lb 0.003 and two renormalization schemes (NDR and HV) are given in 



Table 4: A5 = 1 Wilson coefficients at ^ = rric 
effective flavours. 2/1 = 2/2 = 0. 



1.3 GeV for rrit = 165 GeV and / = 3 





a(1I = 

'"■MS 


290 MeV 


k^ = 

MS 


340 MeV 


^^MS - 


390 MeV 


Scheme 


NDR 


HV 


NDR 


HV 


NDR 


HV 


2/3 


0.027 


0.030 


0.030 


0.034 


0.033 


0.038 


2/4 


-0.054 


-0.056 


-0.059 


-0.061 


-0.064 


-0.067 


2/5 


0.006 


0.015 


0.005 


0.016 


0.003 


0.017 


2/6 


-0.082 


-0.074 


-0.092 


-0.083 


-0.105 


-0.093 


2/7/a 


-0.038 


-0.037 


-0.037 


-0.036 


-0.037 


-0.034 


Vfi/oi 


0.118 


0.127 


0.134 


0.143 


0.152 


0.161 


2/9/a 


-1.410 


-1.410 


-1.437 


-1.437 


-1.466 


-1.466 


2/io/a 


0.496 


0.502 


0.539 


0.546 


0.585 


0.593 



It is customary in phenomenological applications to take Re^o and lo from experiment, 



I.e. 



Re^o = 3.33 • 10"^ GeV, u; = 0.045, (5.11) 

where the last relation reflects the so-called A/ = 1/2 rule. This strategy avoids to a large 
extent the hadronic uncertainties in the real parts of the isospin amplitudes Aj. 

The sum in (|5.7| ) and (5^) runs over all contributing operators. is fully dominated 

by electroweak penguin contributions. on the other hand is governed by QCD penguin 

contributions which are suppressed by isospin breaking in the quark masses (m^ 7^ m-d). The 
latter effect is described by 
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n,„,=i(I^=^. (5.12) 

For ^l-q+T]' we will first set 

Qri+n' = 0.25 , (5.13) 

which is in the ball park of the values obtained in the approach p^] and in chiral 



perturbation theory |65, 67 1. i^r^+r;' is independent of mt- We will investigate the sensitivity 



of e' /e to i7r)+r)' later on. 

5.3 Hadronic Matrix Elements 

The main source of uncertainty in the calculation of e' /e are the hadronic matrix elements 
{Qi)i. They generally depend on the renormalization scale /i and on the scheme used to 
renormalize the operators Qi. These two dependences are canceled by those present in the 
Wilson coefficients so that the resulting physical e' /e does not (in principle) depend on 
^ and on the renormalization scheme of the operators. Unfortunately, the accuracy of the 
present non-perturbative methods used to evalutate {Qi)i is not sufficient to have the ^ and 
scheme dependences of {Qi)i fully under control. We believe that this situation will change 
once the lattice calculations and QCD sum rule calculations improve. A brief review of the 
existing methods including most recent developments will be given below. 

In view of this situation it has been suggested in [^] to determine as many matrix elements 
{Qi)i as possible from the leading CP conserving K vrvr decays, for which the experimental 
data is summarized in (5.11). To this end it turned out to be very convenient to determine 



{Qi)i in the three-flavour effective theory at a scale fi « rric- The details of this approach 
will not be discussed here. It sufficies to say that this method allows to determine only the 
matrix elements of the {V — A) f^i {V — A) operators. For the central value of ImXt these 
operators give a negative contribution to e' /e of about —2.5 • 10~^. This shows that these 
operators are only relevant if e' /e is below 1 • 10~^. Unfortunately the matrix elements of the 
dominant {V — A)<^ {V + A) operators cannot be determined by the CP conserving data and 
one has to use non-perturbative methods to estimate them. 

Concerning the {V — A) ^ {V + A) operators Q5 — Qs, it is customary to express their 

(1/2) (3/2) 

matrix elements {Qi)i in terms of non-perturbative parameters and as follows: 

{Q^)o ^ Bf"^ , ^ Bf"^ . (5.14) 

The label "vac" stands for the vacuum insertion estimate of the hadronic matrix elements in 
question for which B^^^"^"^ = Bf'^'^^ = 1. 
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As the numerical analysis in |74] shows e' je is only weakly sensitive to the values of the 



parameters B- 



(1/2) ^(1/2) fj(l/2) jj(l/2) 



BV''^' and B^^"^^ as long as their absolute values are not 



substantially larger than 1. As in [74 1 our strategy is to set 



B 



(1/2) 
3,7,8 



(mc) = 1, 



B 



(3/2) 



and to treat B^^"^^ (nic 



and B, 



The approach in 
tern 

(3/2) 



MS' 



Bi'/'\m., 



Bf''\m^ 
and B^^'^\mc) as free parameters. 

allows then in a good approximation to express e' /e or equivalently 
Fgi in terms of A^, mt, and the two non-perturbative parameters Bq'^'^^ 

which cannot be fixed by the CP conserving data. 



(5.15) 



5.4 An Analytic Formula for e'/e 



As shown in |^^, it is possible to cast the formal expressions for e'/e in (|5.5| )~(5. 



into an 



analytic formula which exhibits the nit dependence together with the dependence on nis, 
A^, Bq'^'^^ and B^^'^\ To this end the approach for hadronic matrix elements presented 
above is used and ^lr]+rj' is set to 0.25. The analytic formula given below, while being rather 
accurate, exhibits various features which are not transparent in a pure numerical analysis. It 
can be used in phenomenological applications if one is satisfied with a few percent accuracy. 
Needless to say, in the numerical analysis presented below we have used exact expressions. 
In this formulation the function is given simply as follows {xt = m^/M^): 



Fe' 



Po + Px Xoixt) + Py Yoixt) + Pz Zo{xt) + Pe Eoixt). 



with the mt-dependent functions given in subsection 2.4. 

The coefficients Pi are given in terms of Bq^^'^^ 
ms{mc) as follows: 



B. 



(3/2) 



p _ ^(0) , ^(6) p , ^(8) p 



where 



Re, 



B, 



(1/2) 



137 MeV 



1 2 



ms(mc) + md(mc 



Rh = B., 



(3/2) 



137 MeV 



n2 



ms(mc) + md(mc 



(5.16) 

(rric) and 
(5.17) 

(5.18) 



(4^ 
MS' 



The Pi are renormalization scale and scheme independent. They depend, however, on A 
In table | we give the numerical values of rf'\ rp'* and r^^^ for different values of at 
/i = rric in the NDR renormalization scheme |3^]. Actually at NLO only tq coefficients 
are renormalization scheme dependent. The last row gives them in the HV scheme. The 
inspection of table ^ shows that the terms involving Tq^^ and r^'* dominate the ratio e'/e. 
Moreover, the function ZQ{xt) representing a gauge invariant combination of Z^- and 7- 
penguins grows rapidly with rrit and due to < these contributions suppress e'/e strongly 
for large mt IgI, ^ . 
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Table 5: Coefficients in the formula (5.17) for various Ai-i^ in the NDR scheme. The last row 

M V MS 

gives the ro coefficients in the HV scheme. 





Ag^ = 290MeV 




1 = 340MeV 




l = 390MeV 


i 


JO) 

i 


i 


i 


i 


i 


i 


JO) 
i 


J6) 
i 


i 




^0 
Eq 


-2.771 
0.532 
0.396 
0.354 
0.182 


9.779 
0.017 
0.072 
-0.013 
-1.144 


1.429 



-9.404 
0.411 


-2.811 
0.518 
0.381 
0.409 
0.167 


11.127 
0.021 
0.079 
-0.015 
-1.254 


1.267 



-10.230 
0.461 


-2.849 
0.506 
0.367 
0.470 
0.153 


12.691 
0.024 
0.087 
-0.017 
-1.375 


1.081 




-11.164 
0.517 





-2.749 


8.596 


1.050 


-2.788 


9.638 


0.871 


-2.825 


10.813 


0.669 



5.5 The Status of m^, B^^^'^\ Bf''^\ fi^+^z and A^^ 



The present status of these parameters has been recently reviewed in details in |35|. Therefore 
our presentation will be very brief. 

5.5.1 nis 

The present values for ms(2GeV) extracted from lattice calculations and QCD sum rules are 

™.(2GeV)H •'""^^""'^^ (Latfce) [34, H 
(124 ± 22) MeV (QCDS) @ 



The value for QCD sum rules is an average over the results given in QCD sum rules 
also allow to derive lower bounds on the strange quark mass. It was found that generally 
ms(2GeV) ^ 100 MeV |^. If these bounds hold, they would rule out the very low strange 
mass values found in unquenched lattice QCD simulations given above. 

Finally, one should also mention the very recent determination of the strange mass from 
the hadronic r-spectral function [38, ms(2GeV) = (I70I55) MeV. We observe that 



the central value is much larger than the corresponding results given above although the 
error is still large. In the future, however, improved experimental statistics and a better 
understanding of perturbative QCD corrections should make the determination of from 
the r-spectral function competitive to the other methods. On the other hand a very recent 
estimate using new r-like (/)-meson sum rules gives ms(2GeV) = (136 ± 16) MeV []90| . 

We conclude that the error on rris is still rather large. In our numerical analysis of s' /e, 
where m-g is evaluated at the scale mc, we will set 

ms(mc) = (130 ± 25) MeV , (5.20) 
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roughly corresponding to ms{2 GeV) obtained in the lattice approach. 



5.5.2 and 



The values for Bq and obtained in various approaches are collected in table |6|. The 

lattice results have been obtained at /i = 2 GeV. The results in the large-N approach and 
the chiral quark model correspond to scales below 1 GeV. However, as a detailed numerical 
analysis in |74| showed, Bq'^'^^ and B^^"^^ are only weakly dependent on fi. Consequently the 
comparison of these parameters obtained in different approaches at different fi is meaningful. 

Next, the values coming from lattice and chiral quark model are given in the NDR renor- 
malization scheme. The corresponding values in the HV scheme can be found using approx- 
imate relations [pal 



(1/2)^ 



1.2(i?(^/'))NDR, 



(B 



(3/2). 



HV ~ 



(3/2) ^ 



jNDR- 



(5.21) 



The results in the large-N approach are unfortunately not sensitive to the renormalization 
scheme. 



Concerning the lattice results for B^^'^\ the old results read B)^'^q'^\2 GeV) = 1.0 ± 0.2 



(1/2), 



5,6 

94]. More accurate estimates for B^^"^^ have been given in ||^ 



B'i''^\2 GeV) 



0.67 ± 0.04 ± 0.05 (quenched) and b''^''^\2 GeV) = 0.76 ± 0.03 ± 0.05 (/ = 2). However, a 
recent work in shows that lattice calculations of Bq'^'^^ are very uncertain and one has 
to conclude that there are no solid predictions for B^^'^^ from the lattice at present. 



Table 6: Results for Bq'^'^^ and Bo"'^' obtained in various approaches. 



(3/2) 



Method 




r(3/2) 


Lattice[39, 91, 
Large— N|]92|, 
ChQM|81| 


37] 
43 1 


0.72 - 1.10 
1.07-1.58 


0.69 - 1.06 
0.42 - 0.64 
0.75 - 0.79 



We observe that most non-perturbative approaches discussed above found B^ below 
unity. The suppression of -Bg below unity is rather modest (at most 20%) in the lattice 
approaches and in the chiral quark model. In the 1/A^ approach B^^"^^ is rather strongly 
suppressed and can be as low as 0.5. 

Concerning i?g^^^^ the situation is worse. As we stated above there is no solid prediction 
for this parameter in the lattice approach. On the other hand while the average value of 
B^^"^^ in the 1/A^ approach is close to 1.0, the chiral quark model gives in the NDR scheme 
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the value for B^^^"^^ as high as 1.33 ib 0.25. Interestingly both approaches give the ratio 
sJ^/^V^f^^^ in the bah park of 1.7. 

Guided by the results presented above and biased to some extent by the results from the 
large-N approach and lattice calculations, we will use in our numerical analysis below Bq'^'^^ 
and i?g^^^^ in the ranges: 

^ Q ^ Q_3^ ^(3/2) ^ ^ Q 2 (5.22) 

keeping always B^^"^^ > B^^'^\ 
5.5.3 Vlr,+rf and 

The dependence of e' /e on fi^+^z can be studied numerically by using the formula (|5.7| ) or 
incorporated approximately into the analytic formula ( 5.16| ) by simply replacing Sj^/^^ with 
an effective parameter 

(i^r)e.=i?r^^^-^^f|^ (5.23) 

A numerical analysis shows that using (1 — r2^_|_^') overestimates the role of Vtri+'ri' ■ In our 
numerical analysis we have incorporated the uncertainty in 17^+^' by increasing the error in 
S^^^^^ from ±0.2 to ±0.3. 



The last estimates of have been done more than ten years ago [qj]-|67] and it is 

desirable to update these analyses which can be summarized by 

= 0.25 ± 0.08 . (5.24) 

In table ^ we summarize the input parameters used in the numerical analysis of e' /e 

MS 



below. The range for in table |^ corresponds roughly to as(Mz) = 0.119 ± 0.003. 



Table 7: Collection of input parameters. We impose 

^(1/2) > ^(3/2) _ 



Quantity 


Central 


Error 


Reference 


Aa 

MS 


340 MeV 


±50 MeV 


gn 


ms(mc) 


130 MeV 


±25 MeV 


See Text 




1.0 


±0.3 


See Text 


r(3/2) 


0.8 


±0.2 


See Text 
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5.6 Numerical Results for e' je 

In order to make predictions for e' je we need the value of ImAj. This can be obtained from 
the standard analysis of the unitarity triangle as discussed in section 4. 

In what follows we will present two types of numerical analyses of e' je which use the 



methods 1 and 2 discussed already in section 4. This analysis is based on |35]. 

Using the first method we find the probability density distributions for e' je in fig. ^. 
From this distribution we deduce the following results: 

,VeJ("ap-^»-MNDR) 
\ (5.2 . lO- (HV) 

The difference between these two results indicates the left over renormalization scheme de- 
pendence. Since, the resulting probability density distributions for e'/e are very asymmetric 
with very long tails towards large values we quote the medians and the 68%(95%) confidence 
level intervals. This means that 68%(95%) of our data can be found inside the corresponding 
error interval and that 50% of our data has smaller e'/e than our median. 

We observe that negative values of e'/e can be excluded at 95% C.L. For completeness we 
quote the mean and the standard deviation for e'/e: 

I (9,1 ±6,2). 10- (NDR) 
I (6.3 ± 4.8) • 10-4 (HV) 

Using the second method and the parameters in table Q we find : 

1.05 • 10-4 < e' Is < 28.8 • 10"^ (NDR). (5.27) 

and 

0.26 • 10-4 < e'/e < 22.0 • 10-" (HV). (5.28) 



We observe that e'/e is generally lower in the HV scheme if the same values for B^^"^^ 

(3/2) 

and Bg are used in both schemes. Since the present non-perturbative methods do not 
have renormalization scheme dependence fully under control we think that such treatment of 
Sj^/^^ and Sj^/^^ is the proper way of estimating scheme dependences at present. Assuming, 
on the other hand, that the values in (5.22) correspond to the NDR scheme and using the 



relation { f).2l\ }, we find for the HV scheme the range 0.58 • 10 ^ < e'/e < 26.9 • 10 ^ which 
is much closer to the NDR result in ( ^.27| ). This exercise shows that it is very desirable to 
have the scheme dependence under control. 

We observe that the most probable values for e'/e in the NDR scheme are in the ball 
park of 1 • 10-'^. They are lower by roughly 30% in the HV scheme if the same values for 



(Bg^^^\ Bg^^^^ ) are used. On the other hand the ranges in ( 5.27] ) and ( 5.28 ) show that for 
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e'/e [10-4] 



Figure 9: Probability density distributions for e' je in NDR and HV schemes. 
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particular choices of the input parameters, values for e' je as high as (2 — 3) • 10^"^ cannot be 
excluded at present. Let us study this in more detail. 

In table taken from |R5|, we show the values of e'/e in units of 10~^ for specific values 



of B^^I^\Bfl^^ and ms{mc) as calculated in the NDR scheme. The corresponding values 
in the HV scheme are lower as discussed above. The fourth column shows the results for 
central values of all remaining parameters. The comparison of the the fourth and the fifth 
column demonstrates how s' /e is increased when is raised from 340 MeV to 390 MeV. 

' MS 

As stated in ( |5.4| ) e' /e is roughly proportional to Aj^. Finally, in the last column maximal 
values of e'/e are given. To this end we have scanned all parameters relevant for the analysis 
of ImXf within one standard deviation and have chosen the highest value of A^ = 390 MeV. 

^ MS 

Comparison of the last two columns demonstrates the impact of the increase of ImAj from 
its central to its maximal value and of the variation of mt. 

Table ^ gives a good insight in the dependence of e' /e on various parameters which is 
roughly described by (|5.4| ). We observe the following hierarchies: 

• The largest uncertainties reside in nis, Bi'/'^ and e'/e increases universally 
by roughly a factor of 2.3 when ms(mc) is changed from 155 MeV to 105 MeV. The 
increase of B^^'^'^ from 1.0 to 1.3 increases e'/e by (55 ± 10)%, depending on mg and 
Bq . The corresponding changes due to -Bg are approximately (40 it 15)%. 

• The combined uncertainty due to ImAt and nit, present both in ImAt and F^/, is ap- 
proximately ±25%. The uncertainty due to mt alone is only ±5%. 

• The uncertainty due to a|^ is approximately ±16%. 

• The uncertainty due to $1^+^/ is approximately ±12%. 

The large sensitivity of e'/e to nis has been known since the analyses in the eighties. In 



the context of the KTeV result this issue has been analyzed in |98|. It has been found that 



(1/2) (3/2) 

provided 2Bq — B^ < 2 the consistency of the Standard Model with the KTeV result 
requires the 2a bound ms(2 GeV) < 110 MeV. Our analysis is compatible with these findings. 

5.7 Summary 

As we have seen, the estimates of e'/e in the Standard Model are typically below the exper- 
imental data. However, as our scanning analysis shows, for suitably chosen parameters, e'/e 
in the Standard Model can be made consistent with data. However, this happens only if all 
relevant parameters are simultaneously close to their extreme values. This is clearly seen in 
table Moreover, the probability density distributions for e'/e in fig. |^ indicates that values 
of e'/e in the ball park of NA31 and KTeV results are rather improbable. 



39 



Table 8: Values of e' je in units of 10 ^ for specific values of Bq^^'^\ B^^'^'^ and m^{'m(.) and 
other parameters as explained in the text. 
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7.0 
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4.4 


5.2 


6.6 
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9.4 


10.9 


14.1 


1.0 


1.0 


130 


5.5 


6.5 


8.5 






155 


3.3 


4.0 


5.2 



Unfortunately, in view of very large hadronic and substantial parametric uncertainties, it 
is impossible to conclude at present whether new physics contributions are indeed required 
to fit the data. Similarly it is difficult to conclude what is precisely the impact of the e' /e- 



data on the CKM matrix. However, as analyzed in |35| there are indications that the lower 
limit on ImAt is improved. The same applies to the lower limits for the branching ratios for 
Kl Ti^vu and Ki — > T:^e^e~ decays discussed in the following sections. 

It is also clear that the e' /e data puts models in which there are new positive contributions 



to £ and negative contibutions to e' in serious difficulties. In particular as analyzed in |35] 
the two Higgs Doublet Model II |9£] can either be ruled out with improved hadronic matrix 
elements or a powerful lower bound on tan/? can be obtained from e'/e. In the Minimal 
Supersymmetric Standard Model, in addition to charged Higgs exchanges in loop diagrams, 
also charginos contribute. For suitable choice of the supersymmetric parameters, the chargino 
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contribution can enhance e' je with respect to the Standard Model expectations 100 |. Yet, 
generaUy the most conspicuous effect of minimal supersymmetry is a depletion of e' je. The 
situation can be different in more general models in which there are more parameters than 
in the two Higgs doublet model II and in the MSSM, in particular new CP violating phases. 
As an example, in more general supersymmetric models e' je can be made consistent with 



experimental findings |101, 102 1. Unfortunately, in view of the large number of free parameters 
such models are not very predictive. 

The future of e' je in the Standard Model and in its extensions depends on the progress 
in the reduction of parametric and hadronic uncertainties. In any case e' je already played 
a decisive role in establishing direct CP violation in nature and its rather large value gives 
additional strong motivation for searching for this phenomenon in cleaner K decays like 
Ki^ tP'vv and Ki^ TT^e^e~, in B decays, in D decays and elsewhere. We now turn to 
discuss some of these topics. 

6 The Decays tt^utj and Ki^ tt^^utj 

6.1 General Remarks 

We will now move to discuss the semileptonic rare FCNC transitions tt^uv and 

K]^ — > -n^uv. Within the Standard Model these decays are loop-induced semileptonic FCNC 
processes determined only by Z'^-penguin and box diagrams and are governed by the single 
function Xo(x() given in ( |2.15| ). 

A particular and very important virtue of -fC — > nvv is their clean theoretical character. 
This is related to the fact that the low energy hadronic matrix elements required are just the 
matrix elements of quark currents between hadron states, which can be extracted from the 
leading (non-rare) semileptonic decays. Other long-distance contributions are negligibly small 



|103| , 104 1 . As a consequence of these features, the scale ambiguities, inherent to perturbative 
QCD, essentially constitute the only theoretical uncertainties present in the analysis of these 
decays. These theoretical uncertainties have been considerably reduced through the inclusion 



of the next-to-leading QCD corrections |105|-||T0 



The investigation of these low energy rare decay processes in conjunction with their the- 
oretical cleanliness, allows to probe, albeit indirectly, high energy scales of the theory and 
in particular to measure Vtd and ImA^ = ImV^*V(d from tt^uu and Ki^ — > tt^uu re- 

spectively. However, the very fact that these processes are based on higher order electroweak 
effects implies that their branching ratios are expected to be very small and not easy to access 
experimentally. 
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6.2 The Decay K+ tv+uu 
6.2.1 The efTective Hamiltonian 

The effective Hamiitonian for — > vr+z^P can be written as 

7ieff = ^ ^ ."n E (K:Kd^NL + ^t:^td^(2;t))(5-<i)v-A(z^«i^/)y-A. (6.1) 
V2 27r sm^ 9w ^J^^ ^ ^ ^ 

The index /=e, fi, r denotes the lepton flavour. The dependence on the charged lepton mass 
resuhing from the box-graph is negligible for the top contribution. In the charm sector this 
is the case only for the electron and the muon but not for the r-lepton. 
The function X{xt) relevant for the top part is given by 

X{xt) = Xo{xt) + ^Xi{xt)=7jx-Xo{xt), rjx = 0.994, (6.2) 

47r 



with the QCD correction [|To|, |lO|, |lO|] 



X^{xt)=Xl{xt) + 8xt^^^^lnx^. (6.3) 



Here x^ = /i^/M^ with fit = 0{mt) and Xi{xt) is a complicated function given in | 10£ , 



108| , 109 1 . The //j-dependence of the last term in ( |6.3D cancels to the considered order the 



//j-dependence of the leading term Xo(xt(/i)). The leftover /x^-dependence in X{xt) is below 
1%. The factor rjx summarizes the NLO corrections represented by the second term in ( |6.2D . 
With rrit = rnt{mt) the QCD factor r]x is practically independent of rrit and A-^^ and is very 
close to unity. 

The expression corresponding to X{xt) in the charm sector is the function X^^. It 



results from the NLO calculation |107] and is given explicitly in [|109(| . The inclusion of 
NLO corrections reduced considerably the large fic dependence (with fic = 0{mc)) present 
in the leading order expressions for the charm contribution |110|. Varying fic in the range 



1 GeV < /Uc < 3 GeV changes Xnl by roughly 24% after the inclusion of NLO corrections to 
be compared with 56% in the leading order. Further details can be found in |107| , ^. The 
impact of the /Uc uncertainties on the resulting branching ratio Br{K^ ■k^uv) is discussed 
below. 

The numerical values for for A* = "^c and several values of and mdmc) can be 



found in |109f| . The net effect of QCD corrections is to suppress the charm contribution by 
roughly 30%. For our purposes we need only 



^'o(^) = ^ 



2^NL + 3"'^NL 

,(4) 



0.42 ± 0.06 (6.4) 



where the error results from the variation of Ai-^ and mdmr). 

MS 
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6.2.2 Deriving the Branching Ratio 



The relevant hadronic matrix element of the weak current {sd)v-A in ( p.l| ) can be extracted 
with the help of isospin symmetry from the leading decay — > Consequently the 

resulting theoretical expression for the branching fraction Br{K^ tt^uv) can be related 
to the experimentally well known quantity Br{K~^ — > Let us demonstrate this. 

The effective Hamiltonian for the tree level decay — > TT^e'^u is given by 



Using isospin symmetry we have 

{7r+\{sd)v-A\K- 



2{7t'>\{su)v-a\K+)- 



(6.5) 



(6.6) 



Consequently neglecting differences in the phase space of these two decays, due to 7^ m^o 
and TJie 7^ 0, we find 



Br{K+ TT+vi>) 
Br{K-^ 



a 



|Ksp27r2 sin^ Gw ^ 



J2 \v:,V,A + VZVtdX{. 



(6.7) 



-e,fj,,T 



6.2.3 Basic Phenomenology 

Using ( |6.7D and including isospin breaking corrections one finds 



\ A5 



X{xt) 



/ReA, 



V A 



3a'^Br(K+ -n^e+u) c 

= rK+ ^ 2 ■ 4n ^ 

zvr^ sm Ww 



4.11 • 10 



-11 



where we have used 



1 



a = — , sin^ Gw = 0.23, Br{K+ yr^e+z^) = 4.82 • 10" 



(6.8) 
(6.9) 

(6.10) 



Here Aj = V*gVid with Ac being real to a very high accuracy. rx+ = 0.901 summarizes isospin 
breaking corrections in relating — > -k^vu to ■K^e'^v. These isospin breaking 

corrections are due to quark mass effects and electroweak radiative corrections and have been 



calculated in |111| ]. Finally i'o(-'^) is given in ( |6.4D . 

Using the improved Wolfenstein parametrization and the approximate formulae (1.22) 
(1.24) we can next put (|6.8|) into a more transparent form Ijl 



Br{K+ ^ t:+vu) = 4.11 • IQ-^^ A'^X'^{xt)- \{cjfif + (^0 - oY 



where 



2 



(6.11) 
(6.12) 
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The measured value of Br{K^ tt^vv) then determines an elhpse in the {Q,ri) plane 
centered at {qq,^) with 



and having the squared axes 



9l=rl t={"A' (6.14) 



a 



where 



^0 



1 



a ■ Br{K~^ — > ir'^vu) 



(6.15) 



4.11 • 10-11 

Note that ro depends only on the top contribution. The departure of from unity measures 
the relative importance of the internal charm contributions. 

The ellipse defined by rg, qq and a given above intersects with the circle ( |1.32| ). This 
allows to determine q and f/ with 

1 

S = -^ 2 

1 — fj"^ 

and consequently 



- ^o^qI + (1 - ^^){rl - ami) h n=jRl- (6.16) 



Rl = l + Rl-2g, (6.17) 
where ry is assumed to be positive. Given g and fj one can determine Vtd- 

Vtd = AX\l-Q-ifi), \Vtd\=AX'^Rt. (6.18) 

The determination of \Vtd\ and of the unitarity triangle requires the knowledge of Vcb (or 
A) and of l^ufe/V'cfel. Both values are subject to theoretical uncertainties present in the 
existing analyses of tree level decays. Whereas the dependence on iKib/Kbl is rather weak, 
the very strong dependence of Br{K^ vr^z^P) on A or Vcb makes a precise prediction 
for this branching ratio difficult at present. We will return to this below. The dependence 
of Br{K^ t:^vi>) on rrii is also strong. However mt is known already within ±4% and 
consequently the related uncertainty in Br[K^ — > n^vv) is substantialy smaller than the 
corresponding uncertainty due to Vcb- 

6.2.4 Numerical Analysis of —>■ tt^vv 

The uncertainties in the prediction for Br(K^ tt^vv) and in the determination of \Vtd\ 
related to the choice of the renormalization scales nt and ^c in the top part and the charm 
part, respectively have been inestigated in To this end the scales /ic and entering mc{fic) 
and mt{nt), respectively, have been varied in the ranges 1 GeV < /j-c ^ 3GeV and 100 GeV < 
fit < 300 GeV. It has been found that including the full next-to- leading corrections reduces 
the uncertainty in the determination of \Vtd\ from ±14% (LO) to ±4.6% (NLO). The main 
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bulk of this theoretical error stems from the charm sector. In the case of Br{K^ tt^vv), 
the theoretical uncertainty due to ^c,t is reduced from ±22% (LO) to ±7% (NLO). 
Scanning the input parameters of table ^ we find 

Br{K+ TT+uv) = (7.9 ± 3.1) • 10"^^ (6.19) 

where the error comes dominantly from the uncertainties in the CKM parameters. 
It is possible to derive an upper bound on Br[K^ t^'^vu) [ p.09| |: 

2 



(6.20) 



where r^s = ^ymB^/rnB^- This equation translates a lower bound on AMg into an upper 
bound on Br{K^ — > n^uv). This bound is very clean and does not involve theoretical 
hadronic uncertainties except for Vsd- Using 



< 0.2 , A< 0.87 , Po{X) < 0.48 , X{xt) < 1.56 , r.^ < 1.2 (6.21) 



we find 



BriK^ vr+i/p)max = 12.2 • 10"^^ . (6.22) 



This limit could be further strengthened with improved input. However, this bound is strong 
enough to indicate a clear conflict with the Standard Model if Br{K^ tt^uv) should be 
measured at 2 • 10"^''. 

6.2.5 \Vtd\ from tt+vv 

Once Br{K^ — > tt^ui') = Br{K^) is measured, \Vtd\ can be extracted subject to various 
uncertainties: 

IT. I = ^OMscaie ± I ± 0.7^ — ± 0.65———— . (6.23) 

iV^dl \Vcb\ rUc Br(K+) 

Taking a{\Vcb\) = 0.002, a{mc) = 100 MeV and a{Br{K+)) = 10% and adding the errors in 
quadrature we find that \Vtd\ can be determined with an accuracy of ±10%. This number 
is increased to ±11% once the uncertainties due to mt, and |^Mfe|/|Vcf,| are taken into 
account. Clearly this determination can be improved although a determination of \Vtd\ with 
an accuracy better than ±5% seems rather unrealistic. 

6.2.6 Summary and Outlook 

The accuracy of the Standard Model prediction for Br(K^ ■n'^uv) has improved consid- 
erably during the last five years. This progress can be traced back to the improved values of 
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nit and \ Vcb\ and to the inclusion of NLO QCD corrections which considerably reduced the 
scale uncertainties in the charm sector. 

Now, what about the experimental status of this decay ? One of the high-lights of 97 was 



the observation by BNL787 collaboration at Brookhaven [112] of one event consistent with 
the signature expected for this decay. The branching ratio: 

BriK+ ^ TT+ui?) = {A.2tli) • 10"^° (6.24) 

has the central value by a factor of 5 above the Standard Model expectation but in view of 
large errors the result is compatible with the Standard Model. The analysis of additional 
data on ir'^vu present on tape at BNL787 should narrow this range in the near future 

considerably. In view of the clean character of this decay a measurement of its branching 
ratio at the level of 2 • 10^^'^ would signal the presence of physics beyond the Standard Model. 



The Standard Model sensitivity is expected to be reached at AGS around the year 2000 |113]. 
Also Fermilab with the Main Injector could measure this decay |114U . 

6.3 The Decay Ki^ tt^uu 
6.3.1 The effective Hamihonian 

The effective Hamiltonian for Ki^ n^vu is given as follows: 



y/2 27r sin^ 



"^eflf = — ^_ , VtsVtdX{xt){sd)v~A{i'v)v~A + h.c. , (6.25) 



where the function X{xt), present already in — > n'^i'iy, includes NLO corrections and is 
given in (|6.2|) . 

As we will demonstrate shortly, i^L — > tt^i^u proceeds in the Standard Model almost 
entirely through direct CP violation |115U . It is completely dominated by short-distance loop 
diagrams with top quark exchanges. The charm contribution can be fully neglected and the 
theoretical uncertainties present in -k^uu due to mc, fic and A-^ are absent here. 

Consequently the rare decay Ki^ tt^vv is even cleaner than — > -k'^vv and is very well 
suited for the determination of the Wolfenstein parameter r/ and in particular ImAt. 

It is usually stated in the literature that the decay Ki^ -n^vD is dominated by direct 
CP violation. Now the standard definition of the direct CP violation requires the presence of 
strong phases which are completely negligible in i^L ir^vv. Consequently the violation of 
CP symmetry in i^L tt^vv arises through the interference between — mixing and 
the decay amplitude. This type of CP violation is often called mixing-induced CP violation. 



However, as already pointed out by Littenberg [115] and demonstrated explictly in a moment. 



the contribution of CP violation to i^L tt^vv via mixing alone is tiny. It gives 
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Br{Ki^ TT^vv) ~ 2 • 10~^^. Consequently, in this sence, CP violation in Ki^ -k^vv with 
Br[Ki^ — > Ti^uv) = 0(10^^^) is a manifestation of CP violation in the decay and as such 
deserves the name of direct CP violation. In other words the difference in the magnitude of 
CP violation in K]^ — > vrvr (e) and ETl — > tt^vv is a signal of direct CP violation and measuring 
Ki^ — > -K^vi? at the expected level would be another signal of this phenomenon. More details 



on this issue can be found in [116, 117, 118]. 



6.3.2 Deriving the Branching Ratio 

Let us derive the basic formula for Br(Ki^ -k^vv) in a manner analogous to the one for 
Br{K~^ — > TT'^vv). To this end we consider one neutrino flavour and define the complex 
function: 

P = ."n ^tsytdX{xt). (6.26) 

V2 2tt sm^ 9w 

Then the effective Hamiltonian in ( |6.25| ) can be written as 

Wcff = F{sd)v-A{yy)v~A + F*{ds)v^A{yv)v-A ■ (6.27) 
Now, from ( |3.8D we have 

KL = l=[(l + e)K' + (l-e)K'] (6.28) 
where we have neglected | e \^<^ 1. Thus the amplitude for Ki ti^vv is given by 
A{Kl ^ TT^yu) = \F{l+e){A{sd)v-A\K^)+F*{l - E){A{ds)v-A\K^)\ {yy)v-A- 



V2 

Recalling 



(6.29) 

CP\K^) = = (6.30) 



we have 

{A{ds)v-A\K'^) = -{A{sd)v-A\K'^). (6.31) 

where the minus sign is crucial for the subsequent steps. 
Thus we can write 

A{Kl ^ vr^P) = i= [F{1 + e) - F*{1 - e)] (^0|(s(i)y-A|i^°)(^^z^)y-A. (6.32) 

Now the terms e can be safely neglected in comparision with unity, which implies that the 
indirect CP violation (CP violation in the — mixing) is negligible in this decay. We 
have then 

F(l + £-) - F*(l - £-) = -I . Im(y,:y,,) • X{xt). (6.33) 

V2 7rsm fc)w 
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Consequently using isospin relation 

(7r°|(Js)y_A|i^°) = {A{su)v-a\K+) (6.34) 

together with ( |6.5| ) and taking into account the difference in the lifetimes of and we 
have after summation over three neutrino flavours 



a 



Br{K^ 



t{K+) |14sP27r2 sin^ Bw 



[ImXt ■ Xixt)y 



(6.35) 



where \t = V^lVtd 



6.3.3 Master Formulae for Br{Ki^ tt^uv) 

Using (|05|) we can write Br{Ki^ — > -k^uv) simply as follows 

/ Tm A \ ' 



-K+ = 1.80 • 10 



-10 



(6.36) 
(6.37) 



with K-^. given in (|6.9| ) and = 0.944 summarizing isospin breaking corrections in relating 

i^L TT^VV to 



111]. 



Using the Wolfenstein parametrization and ( |6.2| ) we can rewrite ( |6.36 ) as 



Br{Ki^ 



3.0 • 10 



-11 



0.39 



mt{mt) 
170 Gey 



2.3 r 



Vr 



cb 



0.040 



(6.38) 



The determination of r/ using Br{Ki^ -k^vv) requires the knowledge of Vcb and nit. 
The very strong dependence on V^b or A makes a precise prediction for this branching ratio 
difficult at present. 



On the other hand inverting ( |6.36D and using ( |6.2D one finds [118]: 
ImAt = 1.36 • 10 



"170 GeV 


1.15 


'Br{Ki, TT^uu) 


. mt(mt) . 




3 • 10-11 



1/2 



(6.39) 



without any uncertainty in \Vcb\- ( |6.39|) offers the cleanest method to measure ImA^; even 
better than the CP asymmetries in B decays discussed briefly in section 8. 

6.3.4 Numerical Analysis of Kj^ — > -k^vu 

The /Uj-uncertainties present in the function X{xt) have already been discussed in connection 
with — > TT'^i'V'. After the inclusion of NLO corrections they are so small that they can be 
neglected for all practical purposes. Scanning the input parameters of table ^ we find 

Br{Ki, TT^uu) = (2.8 ± 1.1) • 10"^^ (6.40) 

where the error comes dominantly from the uncertainties in the CKM parameters. 
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6.3.5 Summary and Outlook 



The accuracy of the Standard Model prediction for Br{Ki, tt^uu) has improved consid- 
erably during the last five years. This progress can be traced back mainly to the improved 
values of mt and \Vcb\ and to some extent to the inclusion of NLO QCD corrections. 
The present upper bound on Br{Ki^ tt^vt?) from FNAL experiment E799 |119] is 



Br{Ki^ 



•K^vv) < 1.6 • 10"^ 



(6.41) 



This is about five orders of magnitude above the Standard Model expectation ( |6.40D . More- 
over this bound is substantially weaker than the model independent bound Jll6| ] from isospin 
symmetry: 

Br{Ki^ TT^uu) < 4.4 • Br{K+ tt'^vu) (6.42) 
which through ( |6.24 ) gives 

Br{Ki, TT^uv) < 6.1 • 10"^ (6.43) 

Now FNAL-E799 expects to reach the accuracy 0(10"^) and a very interesting new 
experiment at Brookhaven (BNL E926) |113| ] expects to reach the single event sensitivity 
2 • 10~^^ allowing a 10% measurement of the expected branching ratio. There are furthermore 



plans to measure this gold-plated decay with comparable sensitivity at Fermilab |12C] and 



KEK [121]. 




(0,0) 



(1,0) 



Figure 10: Unitarity triangle from K — > vrz^i^. 
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6.4 Unitarity Triangle and sin 2/3 from K txvv 

Tlie measurement of Br{K^ tt^vv) and Br{Ki^ ir^vu) can determine tlie unitarity 
triangle completely, (see fig. |l^), provided mt and Vcb are known |122| ]. Using these two 
branching ratios simultaneously allows to eliminate from the analysis which removes 



a considerable uncertainty. Indeed it is evident from ( |6.8D and ( S.36| ) that, given Br{K'^ 
TT^vv) and Br{K\^ tt^vv), one can extract both ImAf and ReA^. One finds |122| , |2[ 



ImAt = A^^-^ ReAt = -A^^ — \ ' , (6-44) 

X{xt) X{xt) 

where we have defined the "reduced" branching ratios 

4.11-10-11 1.80-10-10 ^ ' 

Using next the expressions for ImAt, ReAt and ReAc given in ( 1.22| )-(1.24) we find 

Po(X)- Va(i?i-i?2) _ ^2 

' = '^ A^) ' ^ = V^A^Xix,) ^'-''^ 

with a defined in ( [6.121 ). An exact treatment of the CKM matrix shows that the formulae 
( 6.46| ) are rather precise | 122 |. 

Using (3) one finds subsequently |122(| 



(^^,^2) = i_i =cot/?, sin 2/3 = -^^ (6.47) 

1 + r J 



with 

^ ^a{Bi-B2)-Po{X) 
rs[Bi,B2) = Vo- 7= (6.48) 



Thus within the approximation of ( |6.46D sin 2/3 is independent of Vet (or A) and mt. 

It should be stressed that sin 2/3 determined this way depends only on two measurable 
branching ratios and on the function Po{X) which is completely calculable in perturbation 
theory. Consequently this determination is free from any hadronic uncertainties and its 
accuracy can be estimated with a high degree of confidence. 



An extensive numerical analysis of the formulae above has been presented in [122, 118 ]. 
Assuming that the branching ratios are known to within ±10% and mt within ±3 GeV one 



finds the results in table ^ |118|. We observe that respectable determinations of all considered 
quantities except for g can be obtained. Of particular interest are the accurate determinations 
of sin 2/3 and of ImAf. The latter quantity as seen in ( 6.3S| ) can be obtained from ETl tt^vv 
alone and does not require knowledge of V^b- The importance of measuring accurately ImAj 
is evident. It plays a central role in the phenomenology of CP violation in K decays and 
is furthermore equivalent to the Jarlskog parameter Jcp ||l^, the invariant measure of CP 
violation in the Standard Model, Jcp = A(l — A^/2)ImAi. 
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Table 9: Illustrative example of the determination of CKM parameters from K tivv. 





cidKftl) = ±0.002 


= ±0.001 


^{\Vtd\) 


±10% 


±9% 


(j{q) 


±0.16 


±0.12 


a{fi) 


±0.04 


±0.03 


cj(sin 213) 


±0.05 


±0.05 


cr(ImAt) 


±5% 


±5% 



The accuracy to which sin 2(3 can be obtained from K — > tivD is, in the example discussed 
above, comparable to the one expected in determining sin 2/3 from CP asymmetries in B 
decays prior to LHC experiments. In this case sin 2/3 is determined best by measuring CP 
violation in J/tpK^. Using the formula for the corresponding time-integrated CP 

asymmetry one finds an interesting connection between rare K decays and B physics |122|] 



^'•f'^"^'' - -acAB.^Jms)l±^ (6.49) 



which must be satisfied in the Standard Model. Here is a B*^ — B^ parameter. We stress 
that except for Po{X) all quantities in (6.49) can be directly measured in experiment and 



that this relationship is essentially independent of rrit and Vcb- Due to very small theoretical 
uncertainties in ( |6.49| ), this relation is particularly suited for tests of CP violation in the 
Standard Model and offers a powerful tool to probe the physics beyond it. 

7 Express Review of Rare K and B Decays 

7.1 The Decays B Xs^vv 

The decays B Xg^dvv are the theoretically cleanest decays in the field of rare -B-decays. 
They are dominated by the same Z'^-penguin and box diagrams involving top quark exchanges 
which we encountered already in the case of and Ki^ n^vi? except for the 

appropriate change of the external quark fiavours. Since the change of external quark fiavours 
has no impact on the rrit dependence, the latter is fully described by the function X{xt) in 
( |6.2| ) which includes the NLO corrections. The charm contribution is fully neglegible here 
and the resulting effective Hamiltonian is very similar to the one for Ki^ — > vr'^z^P given in 
( 6.25| ). For the decay B ^ Xsvv it reads 

Gy ^ ~ 
^eflf = —7=- — r^r7^'^t}ytsX{xt){hs)v~A{i^T^)v-A + h.c. (7.1) 
V 2 27r sm Bw 



51 



with s replaced by d in the case of B ^ Xdvy. 

The theoretical uncertainties related to the renormalization scale dependence are as in 
i^L — > TT^vv and can be essentially neglected. The same applies to long distance contributions 



considered in [123]. The calculation of the branching fractions for B — > Xg^^uv can be done in 
the spectator model corrected for short distance QCD effects. Normalizing to Br{B Xcev) 
and summing over three neutrino flavours one finds 

Br{B ^ Xsvv) ^ \Vts? X'^jxt) 

Br{B^Xceiy) 47r2 sin^ Gw l^cbP f{z) k{z) ' 
Here f{z) is the phase-space factor for B X^ev with z = m^/m^ and k,{z) = 0.88 [ 124 , 



|125([ is the corresponding QCD correction. The factor k(0) = 0.83 represents the QCD 
correction to the matrix element of the h — > svv transition due to virtual and bremsstrahlung 
contributions. In the case of -B ^ X^uU one has to replace Vts by Vtd which results in a 
decrease of the branching ratio by roughly an order of magnitude. 

Setting Br{B Xccu) = 10.4%, f{z) = 0.54, k{z) = 0.88 and using the values in ( |6loD 
we have 

(7.3) 



5r(i?^X..P) = 3.7.10-^"^*^"' r-tM^^-^° 



\Vcb\^ Ll70GeV 

Taking next, f{z) = 0.54 ± 0.04 and Br{B XcCv) = (10.4 ± 0.4)% and scanning the input 
parameters of table |2| we find 

Br{B Xgvv) = (3.5 ± 0.7) • 10"^ (7.4) 

to be compared with the experimental upper bound: 

Br{B Xsvv) < 7.7 • 10"^ (90% C.L.) (7.5) 



obtained for the first time by ALEPH [12£]. This is only a factor of 20 above the Standard 
Model expectation. Even if the actual measurement of this decay is extremly difficult, all 
efforts should be made to measure it. One should also make attempts to measure Br{B — > 
Xdi^i^). Indeed 

Br{B^Xdi^=\3£ (7n 
Br{B^XgVv) \Vts\^ 

offers the cleanest direct determination of [Vtd[/[Vjs[ as all uncertainties related to mt, f{z) 

and Br[B Xcev) cancel out. 

7.2 The Decays B^^d I'^l' 

The decays B^ d ~^ /^/^ are after B Xg^d^^U the theoretically cleanest decays in the field 
of rare i?-decays. They are dominated by the Z'^-penguin and box diagrams involving top 
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quark exchanges which we encountered already in the case of i? — > Xs^d^i^ except that due to 
charged leptons in the final state the charge flow in the internal lepton line present in the box 
diagram is reversed. This results in a different mt dependence summarized by the function 
Y{xt), the NLO generalization |l06|, |l08|, |lO|] of the function Yo{xt) given in (|2l^ ). The 
charm contributions are fully negligible here and the resulting effective Hamiltonian is given 
for Bg l^l^ as follows: 



a 



V2 2tt sin^ Bw 

with s replaced by d in the case of l^l 
The function Y{x) is given by 



V*^VtsY{xt){hs)v-A{ll)v-A + h.c. 



(7.7) 



Y{xt) = Yo{xt) + -^Yi{xt) = VvYoixt), 
47r 



VY 



1.012 



(7.8) 



where Yi{xt) can be found in [ 106| , 108 , 109]. The leftover ^u^-dependence in Y{xt) is tiny and 
amounts to an uncertainty of ±1% at the level of the branching ratio. With mt = mt{mt) the 
QCD factor rjy depends only very weakly on mt. The dependence on A.-^ can be neglected. 
The branching ratio for Bg — > is given by Jl06|] 



BriBg 



1+1-) = t{Bs) 



r(Bg)^ 

IT 



a 



Att sin 0w 



\ 



mf 



V:M^Y\xt) (7.9) 



ml 



where Bs denotes the flavour eigenstate (bs) and Fb^ is the corresponding decay constant. 
Using ( S.lOl ) and {^^) we find in the case of Bg 



Br{Bs /U+/Z") = 3.5 • 10"^ 



riBg) 



1.6ps 



Wt 



-.2 r 



ts 



0.040 



170 GeV 



3.12 



(7.10) 



_210MeV. 

The main uncertainty in this branching ratio results from the uncertainty in Fb^ ■ Scanning 
the input parameters of table |2| together with t{Bs) = 1.6 ps and Fb^ = (210 it 30) MeV we 
find 

Br{Bs /x+/i") = (3.2 ± 1.5) • 10^^ . (7.11) 

For Bd A*"*"/^" a similar formula holds with obvious replacements of labels (s — > d). Pro- 
vided the decay constants Fb^ and Fb^ will have been calculated reliably by non-perturbative 
methods or measured in leading leptonic decays one day, the rare processes Bg — > fJ-^fJ-^ and 
Bd — > fJ'~^fJ'~ should offer clean determinations of \Vtg\ and \Vtd\- In particular the ratio 



T{Bd) mB,Fl\Vtd\' 



BrjBd ^ /i+/i-) 
Br{Bg ^ /i+/x-) T{Bg) ruB^ F| \Vtg\ 



(7.12) 



having smaller theoretical uncertainties than the separate branching ratios should offer a 
useful measurement of |Vtd|/|Vts|. Since Br{Bd — > iJ-^/J'^) = 0(10^^'^) this is, however, a very 
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difficult task. For Bg — > r^r^ and Bs e^e~ one expects branching ratios 0(10^^) and 
0(10"^^), respectively, with the corresponding branching ratios for S^^-decays by one order 
of magnitude smaller. 

The bounds on Bg ^ — > II are still many orders of magnitude away from Standard Model 
expectations. The best bounds come from CDF pi. One has: 



Br{Bs ^+11-) < 2.6 • 10-6 (95%C7.L.) (7.13) 

and Br{Bd /U^/U^) < 8.6 • 10^''. CDF should reach in Run II the sensitivity of 1 • 10~^ 
and 4 • 10^*^ for B^ fj.fi and Bs — > fj.p-, respectively. It is hoped that these decays will be 
observed at LHC-B. The experimental status B ^ t^t^ and its usefulness in tests of the 



physics beyond the Standard Model is discussed in |125]. 



7.3 B X,7 and B X^l^l' 

In view of space limitations I will be very brief on these two decays. 

A lot of efforts have been put into predicting the branching ratio for the inclusive radiative 
decay B Xsj including NLO QCD corrections and higher order electroweak corrections. 
The relevant references can be found in |||, ^l|, where also theoretical details are given. The 
final result of these efforts can be summarized by 

Br{B Xs-f)th = (3.30 ± 0.15(scale) ± 0.26(par)) • 10"^ (7.14) 

where the first error represents residual scale dependences and the second error is due to 
uncertainties in input parameters. The main achievement is the reduction of the scale depen- 
dence through NLO calculations, in particular those given in p9t| and pH]. In the leading 
order the corresponding error would be roughly ±0.6. 



The theoretical result in ( [7.14 ) should be compared with experimental data: 



, (3.15 ± 0.35 ± 0.41) • lO-'' , CLEO 
Br{B ^ X,7 exp = <^ / 7.15 

' (3.11 ±0.80 ±0.72) • lO"'' , ALEPH, 

which implies the combined branching ratio: 

Br{B ^ X,7)exp = (3.14 ± 0.48) • 10"^ . (7.16) 

Clearly, the Standard Model result agrees well with the data. In order to see whether any 
new physics can be seen in this decay, the theoretical and in particular experimental errors 
should be reduced. This is certainly a very difficult task. 

The rare decays B have been the subject of many theoretical studies. It is 

clear that once these decays have been observed, they will offer useful tests of the Standard 



Model and of its extentions. Most recent reviews can be found in 131 , 132 1 
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7.4 Kl 7r°e+e- 

There are three contributions to this decay: CP conserving, indirectly CP violating and 
directly CP violating. Unfortunately out of these three contributions only the directly CP 



violating can be calculated reliably. Including NLO corrections 133 | and scanning the input 
parameters of table |2| we find 

Br{Ki^ 7r°e+e-)dir = (4.6 ± 1.8) • lO'^^ , (7.17) 

where the errors come dominantly from the uncertainties in the CKM parameters. The 



remaining two contributions to this decay are plagued by theoretical uncertainties |134]. 
They are expected to be 0(10^^^) but generally smaller than Br{Ki^ — > 7r'^e^e~)dir- This 
implies that within the Standard Model Br{Ki^ — > Tr^e^e~) is expected to be at most 10^^^. 
Experimentally we have the bound |135|] 



Br{Ki^ 7r°e+e-) < 4.3 • 10"^ (7.18) 
and considerable improvements are expected in the coming years. 

7.5 

The i^L Ai^Ai^ branching ratio can be decomposed generally as follows: 

BR{Ki^ fi+fi-) = \ReAf + \lmAf , (7.19) 

where HeA denotes the dispersive contribution and 1mA the absorptive one. The latter 
contribution can be determined in a model independent way from the Kl — > 77 branching 
ratio. The resulting |Im^p is very close to the experimental branching ratio Br{Ki^ — > 
fj-^fj-^) = (7.2 lb 0.5) • 10^^ |136| ] so that |Re^p is substantially smaller and extracted to be 

m 

iReAxpP < 5.6 • 10"^° (90% C.L.). (7.20) 
Now ReA can be decomposed as 

Re^ = Re^LD + Re^sD , (7.21) 

with 

\ReAsB\^ = BriKi^^ fi+fi-)sB (7.22) 
representing the short-distance contribution which can be calculated reliably. An improved 



estimate of the long-distance contribution Re^Lo has been recently presented in | 137 ] 

iRe^LDl < 2.9 • 10"^ (90% C.L.). (7.23) 
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Together with ( |7.20| ) this gives 



Br{KL ^ /U+^-)sD < 2.8 • 10" 



(7.24) 



This result is very close to the one presented by Gomez Dumm and Pich |13^ ]. More pesimistic 



view on the extraction of the short distance part from Br{Ki^ i^i^fj, ) can be found in [139]. 



The bound in ( 7.24 ) should be compared with the short distance contribution within the 
Standard Model for which we find 



Sr(i^L ^ /^i+/U")sD = (8.7 ± 3.6) • 10 



-10 



(7.25) 



This implies that there is a considerable room for new physics contributions. We will return 
to this point in section 9. Reviews of rare K decays are hsted in [pO| . 



8 Express Review of CP Violation in B Decays 
8.1 CP- Asymmetries in B-Decays: General Picture 

CP violation in B-decays is certainly one of the most important targets of B-factories and 
of dedicated B-experiments at hadron facilities. It is well known that CP violating effects 
are expected to occur in a large number of channels at a level attainable at forthcoming 
experiments. Moreover there exist channels which offer the determination of CKM phases 
essentially without any hadronic uncertainties. Since extensive reviews on CP violation in 
B decays can be found in the literature |141, 142, ^ and I am running out of space, let me 



concentrate only on a few points beginning with a quick review of classic methods for the 
determination of the angles a, (3 and 7 in the unitarity triangle. 

The classic determination of a by means of the time dependent CP asymmetry in the 
decay B'^ — > 7r~^7:~ is affected by the "QCD penguin pollution" which has to be taken care 
of in order to extract a. The recent CLEG results for penguin dominated decays indicate 
that this pollution could be substantial as stressed in particular in |143|. The most popular 



strategy to deal with this "penguin problem" is the isospin analysis of Gronau and London 
1 144]. It requires however the measurement of Br{B^ vr'^vr'^) which is expected to be below 



10 ^: a very difficult experimental task. For this reason several, rather involved, strategies 



]145] have been proposed which avoid the use oi B^ tt^tt^ in conjunction with acp(7r+7r^, t). 
They are reviewed in ]|3| . It is to be seen which of these methods will eventually allow us to 
measure a with a respectable precision. It is however clear that the determination of this 
angle is a real challenge for both theorists and experimentalists. 

The CP-asymmetry in the decay B^ ipKs allows in the Standard Model a direct 
measurement of the angle /? in the unitarity triangle without any theoretical uncertainties 
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Ill4q] . Of considerable interest |142| , |147] ] is also the pure penguin decay 4>Ks, which is 

expected to be sensitive to physics beyond the Standard Model. Comparision of /3 extracted 
from Bd — > (pKs with the one from B^ — > ijjKs should be important in this respect. An 



analogue of B^i — > 'ipKs in S^-decays is Bg — > tpcj). The CP asymmetry measures here rj | 14S | 
in the Wolfenstein parametrization. It is very small, however, and this fact makes it a good 
place to look for the physics beyond the Standard Model. In particular the CP violation 
in B^ — Bg mixing from new sources beyond the Standard Model should be probed in this 
decay. 

The two theoretically cleanest methods for the determination of 7 are: i) the full time 



dependent analysis of Bg — > DfK and Bg Dg [149] and ii) the well known triangle 



construction due to Gronau and Wyler [150| which uses six decay rates B^ — > D^^pK^, 



B+ D^K+, D^K+ and B D^K', D^K'. Both methods are unaffected by penguin 
contributions. The first method is experimentally very challenging because of the expected 
large B^ — B^ mixing. The second method is problematic because of the small branching 
ratios of the colour supressed channel B~^ — > D^K~^ and its charge conjugate, giving a rather 
squashed triangle and thereby making the extraction of 7 very difficult. Variants of the latter 



method which could be more promising have been proposed in [151, 152 1. It appears that 
these methods will give useful results at later stages of CP-B investigations. In particular the 
first method will be feasible only at LHC-B. Other recent strategies for 7 will be mentioned 
below. 

8.2 B°-Decays to CP Eigenstates 

Let us demonstrate some of the statements made above explicitly. 

A time dependent asymmetry in the decay B^ — > / with / being a CP eigenstate is given 

by 

acp{t, f) = A'^UB f) cos(AMt) + ^ /) sin(AMt) (8.1) 

where we have separated the direct CP-violating contributions from those describing mixing- 
induced CP violation: 

In (|8.lD , AM denotes the mass splitting of the physical B^-B^ -Tailing eigenstates. The 
quantity containing essentially all the information needed to evaluate the asymmetries 
(|8.2|) is given by 

e/ = exp(i20M)^||^yy (8.3) 
with (j)M denoting the weak phase in the B — B mixing and A{B — > /) the decay amplitude. 
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Generally several decay mechanisms with different weak and strong phases can contribute 
to A{B /). These are tree diagram (current-current) contributions, QCD penguin con- 
tributions and electroweak penguin contributions. If they contribute with similar strength 
to a given decay amplitude the resulting CP asymmetries suffer from hadronic uncertainies 
related to matrix elements of the relevant operators Qi. 

An interesting case arises when a single mechanism dominates the decay amplitude or the 
contributing mechanisms have the same weak phases. Then 

= exp(i2</>M) exp(-i2(^D), I P= 1 (8.4) 

where (^_d is the weak phase in the decay amplitude. In this particular case the hadronic 
matrix elements drop out, the direct CP violating contribution vanishes and the mixing- 
induced CP asymmetry is given entirely in terms of the weak phases (j)M and (pD- In particular 
the time integrated asymmetry is given by 

acp{f) = ±sm{24>D - 20m) (8.5) 

where it refers to / being a CP = it eigenstate and x^^s are the B'^ s~Bd s mixing parameters. 
If a single tree diagram dominates, the factor sm{2(p£) — 2(j)M) can be calculated by using 

^ / ^ ^ / (R(^^ 

4>D = { , 0A/ = < „ (8.6) 

[ [ 5° 

where we have indicated the basic transition of the b-quark into a lighter quark. On the 
other hand if the penguin diagram with internal top exchange dominates one has 





D={ ' . <t>M = { J (8.7) 



These rules have been obtained using Wolfenstein parametrization in the leading order. Let 
us practice these formulae. Assuming that Bd ipKs and B^ 7r"^7r~ are dominated by 
tree diagrams with b ^ c and b ^ u transitions respectively we readly find 

acpms) = -sm{2P)-^, (8.8) 

acp{7T+TT~) = -sm{2a)-^. (8.9) 

l + xj 

Now in the case of B^ ipKs the penguin diagrams have to a very good approximation 
the same phase {(pD =0) as the tree contribution and moreover are Zweig suppressed. Con- 
sequently (|8.8|) is very accurate. This is not the case for Bd tt^tt^ where the penguin 
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contribution could be substantial. Heaving weak phase (1>d = —P, which differs from the tree 
phase cjyD = 7, this penguin contribution changes effectively ( |8.9| ) to 



acpiTT+TT") = -sm{2a + 9p)—^ (8.10) 



where is a function of (3 and hadronic parameters. The isospin analysis [|144| mentioned 
before is supposed to determine 6p so that a can be extracted from acp (7r''~7r~). 

Similarly the pure penguin dominated decay (pKs is governed by the b ^ s penguin 

with internal top exchange which implies that in this decay the angle /3 is measured. The 
accuracy of this measurement is a bit lower than using — > ipKs as penguins with internal 
u and c exchanges may introduce a small pollution. 

Finally we can consider the asymmetry in Bg — > 4^<j3, an analog of Bj, ipKg. In the 
leading order of the Wolfenstein parametrization the asymmetry acp{il^(j)) vanishes. Including 



higher order terms in A one finds |145 



acpW) = 2A%— ^ (8.11) 
I + xi 



where A and r] are the Wolfenstein parameters. 



8.3 Recent Developments 



All this has been known already for some time and is well documented in the literature. The 
most recent developments are related to the extraction of the angle 7 from the decays B — > PP 



(P=pseudoscalar) and their charge conjugates | 153 |- |15q] . Some of these modes have been 
observed by the CLEO collaboration ]157| ]. In the future they should allow us to obtain 
direct information on 7 at S-factories (BaBar, BELLE, CLEO III) (for interesting feasibility 



studies, see [154, |155| , |131[| ). At present, there are only experimental results available for the 



combined branching ratios of these modes, i.e. averaged over decay and its charge conjugate, 
suffering from large hadronic uncertainties. 

There has been large activity in this field during the last two years. The main issues 
here are the final state interactions, SU(3) symmetry breaking effects and the importance 
of electroweak penguin contributions. Several interesting ideas have been put forward to 



extract the angle 7 in spite of large hadronic uncertainties in B ^ ttK decays |153|, 154] 



Also other B PP decays have been investigated. As this field became rather technical, I 



decided not to include it in these lectures. A subset of relevant papers is listed in [153, 154, 



156, |158| , 159, 16C], where further references can be found. In particular in [156, [L59[ [ general 
parametrizations for the study of the final state interactions, SU(3) symmetry breaking effects 
and the importance of electroweak penguin contributions have been presented. Moreover, 
upper bounds on the latter contributions following from SU(3) symmetry have been derived 
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|16C]. Recent reviews can be found in [161] and 162 |. New strategies for 7 which include 
Bg ipKs ^-ncl Bs — > K~ have been suggested very recently in |163f| . 

There is no doubt that these new ideas will be helpful in the future. They are, however, 
rather demanding for experimentalist as often several branching ratios have to be studied 
simultaneously and each has to be measured precisely in order to obtain an acceptable mea- 
surement of 7. On the other hand various suggested bounds on 7 may either exclude the 



region around 90° |15S] or give an improved lower bound on it 160 , 162 , 164] which would 



remove a large portion of the allowed range from the analysis of the unitarity triangle. In this 



context it has been pointed out in 1165] (see also [164]) that generally charmless hadronic B 
decay results from CLEO seem to prefer negative values of cos 7 which is not the case in the 
standard analysis of section 4. 

Finally I would like to mention a recent interesting paper of Lenz, Nierste and Ostermaier 
|166| ], where inclusive direct CP-asymmetries in charmless S^-decays including QCD effects 
have been studied. These asymmetries should offer additional useful means to constrain the 
unitarity triangle. 

8.4 CP- Asymmetries in i?-Decays versus K — > -nvu 

Let us next compare the potentials of the CP asymmetries in determining the parameters of 
the Standard Model with those of the cleanest rare /C-decays: Ki^ tt^vD and -k'^uu. 
Measuring sin 2a and sin 2/3 from CP asymmetries in B decays allows, in principle, to fix 



the parameters r] and which can be expressed as [167] 



r_ (sin 2a) -|- r_|_ (sin 2/3) 
l + r\{sm 2/3) 



^=^-^^7-^7^5^' ^= l-f?r+(sin2/3), (8.12) 



where r±{z) = (1 it Vl — z^)l z. In general the calculation of q and r/ from sin 2a and sin 2/3 
involves discrete ambiguities. As described in |167] ] they can be resolved by using further 



information, e.g. bounds on iKifc/Fcfc], so that eventually the solution (|8.12D is singled out. 

Let us then consider two scenarios of the measurements of CP asymmetries in B^ vr+vr^ 
and — > J/ij^Ko,, expressed in terms of sin 2a and sin 2/3: 

sin 2a = 0.40 ±0.10, sin 2/3 = 0.70 ± 0.06 (scenario I) (8.13) 

sin 2a = 0.40 ± 0.04 , sin 2(3 = 0.70 ± 0.02 (scenario II) . (8.14) 

Scenario I corresponds to the accuracy being aimed for at i3-factories and HERA-B prior to 
the LHC era. An improved precision can be anticipated from LHC experiments, which we 
illustrate with the scenario II. We assume that the problems with the determination of a will 
be solved somehow. 



60 



Table 10: Illustrative example of the determination of CKM parameters from K tivv and 
B-decays. We use a{\Vcb\) = ±0.002(0.001). 





K — > TTW 


Scenario I 


Scenario II 


^{\ytd\) 


±10%(9%) 


±5.5%(3.5%) 


±5.0%(2.5%) 


a{Q) 


±0.16(0.12) 


±0.03 


±0.01 


cr(?7) 


±0.04(0.03) 


±0.04 


±0.01 


cj(sin2/3) 


±0.05 


±0.06 


±0.02 


cr(ImAt) 


±5% 


±14%(11%) 


±10%(6%) 



In table IC this way of the determination of the Standard Model parameters is compared 
1 115] with the analogous analysis using Ki^ -k^vu and ■k'^vv which has been pre- 

sented in section 6. As can be seen in table |l^, the CKM determination using K — > Trvv is 
competitive with the one based on CP violation in B decays in scenario I, except for g which 
is less constrained by the rare kaon processes. On the other hand as advertised previously 
ImAf is better determined va. K ^ t^uv even if scenario II is considered. The virtue of the 
comparision of the determinations of various parameters using CP-B asymmetries with the 
determinations in very clean decays K tivv is that any substantial deviations from these 
two determinations would signal new physics beyond the Standard Model. Formula ( |6.49D 
is an example of such a comparison. There are other strategies for determination of the 
unitarity triangle using combinations of CP asymmetries and rare decays. They are reviewed 
in m. 



9 A Brief Look Beyond the Standard Model 

9.1 General Remarks 

We begin the discussion of the Physics beyond the Standard Model with a few general re- 
marks. As the new particles in the extensions of the Standard Model are generally substan- 
tally heavier than W^, the impact of new physics on charged current tree level decays should 
be marginal. On the other hand these new contributions could have in principle an important 
impact on loop induced decays. From these two observations we conclude: 

• New physics should have only marginal impact on the determination of \Vus\, and 

\Vub\. 
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• There is no impact on the calculations of the low energy non-perturbative parameters 
Bi except that new physics can bring new local operators implying new parameters Bi . 

• New physics could have substantial impact on rare and CP violating decays and con- 
sequently on the determination of the unitarity triangle. 

9.2 Classification of New Physics 

Let us then group the extensions of the Standard Model in three classes. 
Class A 

• There are no new complex phases and quark mixing is described by the CKM matrix. 

• There are new contributions to rare and CP violating decays through diagrams involving 
new internal particles. 

These new contributions will have impact on the determination of a, /3, 7, \Vtd\ and At and 
will be signaled by 

• Inconsistencies in the determination of r/) through e, B^ ^ — B^ ^ mixing and rare 
decays. 

• Disagreement of {g, fj) extracted from loop induced decays with {g, fj) extracted using 
CP asymmetries. 

Examples are two Higgs doublet model II and the constrained MSSM. 
Class B 

• Quark mixing is described by the CKM matrix. 

• There are new phases in the new contributions to rare and CP violating decays. 

This kind of new physics will also be signaled by inconsistencies in the (g, f)) plane. However, 
new complication arises. Because of new phases CP violating asymmetries measure generally 
different quantities than a, [3 and 7. For instance the CP asymmetry m. B ^ ipKs will no 
longer measure /3 but 13 + 9np where 6'ivp is a new phase. Strategies for dealling with such 



situation have been developed. See for instance [|116| , 165] and references therein. 

Examples are multi-Higgs models with complex phases in the Higgs sector, general SUSY 
models, models with spontaneous CP violation and left-right symmetric models. 

Class C 

• The unitarity of the three generation CKM matrix does not hold. 
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Examples are four generation models and models with tree level FCNC transitions. If this 
type of physics is present, the unitarity triangle will not close or some inconsistencies in the 
{g, ff) plane take place. 

Clearly in order to sort out which type of new physics is responsible for deviations from 
the Standard Model expectations one has to study many loop induced decays and many CP 



asymmeteries. Some ideas in this direction can be found in |165, llC]. 



9.3 Upper Bounds on K -nvv and Ki^ n^e^e from e'/e and Kl /x+yU 

We have seen in previous sections that the rare kaon decays Ki^ — > -k^vv, -n'^vD 
and ETl ~^ n^e'^e" are governed by Z-penguin diagrams. Within the Standard Model the 
branching ratios for these decays have been found to be 

Br{Ki^^ -K^uv) = (2.8± 1.1) • 10"^S (9.1) 
Br{K+ ^ TT+iyj?) = (7.9±3.1) • lO^^S (9.2) 
-Br(A"L ^ vr°e+e")dir = (4.6 ± 1.8) • 10"^^^ (9.3) 

where the errors come dominantly from the uncertainties in the CKM parameters. The 
branching ratio in ( |9.3D represents the so-called direct CP-violating contribution to Ki^ — > 
■K^e~^e~. The remaining two contributions to this decay, the CP-conserving one and the 



indirect CP-violating one are plagued by theoretical uncertainties |134]. They are expected 
to be 0(10"^^) but generally smaller than Br{Ki^ — > 7r''e'''e~)dir- This implies that within 
the Standard Model Br{Ki^ — > vr'^e+e") is expected to be at most 10^^-*^. 

In this context a very interesting claim has been made by Colangelo and Isidori |169| ], 
who analyzing rare kaon decays in supersymmetric theories pointed out a possible large 
enhancement of the effective sdZ vertex leading to an enhancement of Br{K^ tt^vv) 
by one order of magnitude and of Br{Ki tt^vv) and Br{Ki^ -K^e^e^) by two orders 
of magnitude relative to the Standard Model expectations. Not surprisingly these results 
brought a lot of excitement among experimentalists. 

Whether substantial enhancements of the branching ratios in question are indeed possible 
in supersymmetric theories is being investigated at present. On the other hand it can be 
shown jl70| ] that in models in which the dominant new effect is an enhanced sdZ vertex, 
enhancements of Br{Ki^ -K^ui?) and Br{Ki^ vr^e^e^) as large as claimed in |169|| are 
already excluded by the existing data on e'/e in spite of large theoretical uncertainties. 
Similarly the large enhancement of Br[K^ t^^vU) can be excluded by the data on e'/e 
and in particular by the present information on the short distance contribution to Ki^ — > 
■ The latter can be bounded by analysing the data on Br(Ki^ fi^fi^) in conjunction 
with improved estimates of long distance dispersive contributions [|137| , |138| . In |16S|] only 
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constraints from Ki^ fi'^ fj,~ , the K^-Ks mass difference AMk and e have been taken into 
account. As e'/e depends sensitively on the size of Z-penguin contributions and generally on 
the size of the effective sdZ vertex it is clear that the inclusion of the constraints from e' /e 
should have an important impact on the bounds for the rare decays in question. I will only 
describe the basic idea of |170| and give numerical results. The relevant expressions can be 



found in this paper. Here we go. 

In the Standard Model Z-penguins are represented by the function Co which enters the 
functions Xq, Yq and Zq. In order to study the effect of an enhanced sdZ vertex one simply 
makes the following replacement in the formulae for e'/e, Ki^ — > /x^/i^ and rare decays in 
question: 

XtCoixt) =^ Zds (9.4) 

where Zds denotes an effective sdZ vertex. The remaining contributions to e' /e, Ki^ /^^^^ 
and rare K decays are evaluated in the Standard model as we assume that they are only 
marginally affected by new physics. We will, however, consider three scenarios for Xt, which 
enters these remaining contributions. 

Indeed there is the possibility that the value of Aj is modified by new contributions to e 



and i?^^ — B'^g mixings. We consider therefore three scenarios: 



• Scenario A: At is taken from the standard analysis of the unitarity triangle 

• Scenario B: ImAt = and ReAt is varied in the full range consistent with the unitarity 
of the CKM matrix. In this scenario CP violation comes entirely from new physics 
contributions. 

• Scenario C: At is varied in the full range consistent with the unitarity of the CKM 
matrix. This means in particular that ImAt can be negative. 



Table 11: Upper bounds for the rare decays Ki^ tt^vv, K]^ — > ir^e^e and ■k'^vv, 
obtained in various scenarios by imposing e' /e > 2.5 • 10~^, in the case ImZds > 0. 



Scenario 


A 


B 


C 


SM 


Br{KL 7r°i/i7)[10-i°] 


0.5 




0.7 


0.4 


Br{KL 7r°e+e-) [10-11] 


0.8 




1.0 


0.7 


Br{K+ 7r+i/p)[10-i°] 


1.8 




2.2 


1.1 



Now Zds is a complex number. ImZds can be best bounded by e'/e. This implies bounds 
for Br{Ki^ — > -n^vU) and Br{Ki^ Tr^e~^e~) which are sensitive functions of ImZds- KeZ^s 
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Table 12: Upper bounds for the rare decays ETl 'k^vv, Ki^ ir^e^e and — > ■k'^uv^ 
obtained in various scenarios by imposing e' /e > 1.5 • 10^^, in the case ImZds > 0. 



Scenario 


A 


B 


C 


SM 


Br{KL 7r°i/i?)[10-^°] 


1.1 




1.2 


0.4 


Br{KL vrOe+e") [10-11] 


1.5 




1.8 


0.7 


Br{K+ 7r+i/p)[10-i°] 


1.9 




2.3 


1.1 



can be bounded by the present information on the short distance contribution to Ki^ l^^f^ - 
This bound imphes a bound on Br{K^ tt^uu). Since Br{K^ tt^uD) depends on both 
ReZ^<. and IvaZ^g also the bound on IvuZ^s from e' /e matters in cases where ImZf^g is very 
enhanced over the Standard Model value. 

The branching ratios Br{K\^ tt^vv) and Br{Ki — > TT^e'^e~) are dominated by (ImZ^^)^. 
Yet, the outcome of this analysis depends sensitively on the sign of ImZg^. Indeed, ImZgd > 
results in the suppression of e' /e and as in the Standard Model the value for e'/e is generally 
below the data substantial enhancements of ImZgd with ImZgd > are not possible. The 
situation changes if new physics reverses the sign of ImZgd so that it becomes negative. Then 
the upper bound on ImZgd is governed by the upper bound on e' /e and with suitable choice 
of hadronic parameters and ImAt (in particular in scenario C) large enhancements of —ImZgd 
and of rare decay branching ratios are possible. The largest branching ratios are found when 
the neutral meson mixing is dominated by new physics contributions which force ImAf to 
be as negative as possible within the unitarity of the CKM matrix. This possibility is quite 
remote. However, if this situation could be realized in some exotic model, then the branching 
ratios in question could be very high. 



In table 11 we show the upper bounds on rare decays for ImZg^ > for three scenarios in 
question and e' /s > 2.5 • 10"^. In table 12 the corresponding bounds for e'/e > 1.5 • 10"^ are 
given. To this end all parameters relevant for e' /e have been scanned in the ranges used in 
section 5. In tables |l3| and |14| the case ImZ^^ < for e'/e < 2.0 • 10'^ and e'/e < 3.0 • 10"^ 
is considered respectively. In the last column we always give the upper bounds obtained in 
the Standard Model. Evidently for positive ImZgd the enhancement of branching ratios are 
moderate but they can be very large when ImZsd < 0. 

Other recent extensive analyses of supersymmetry effects in ^ vrz^P have been presented 



in [116, [171| , 172| where further references can be found. Model independent studies of these 
decays can be found in [116, 172|. The corresponding analyses in various no-supersymmetric 



extensions of the Standard Model are listed in [173|. In particular, enhancement of Br{Ki 
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Table 13: Upper bounds for the rare decays ETl — > 'k^vv, Ki^ ir^e^e and ■k'^uv^ 
obtained in various scenarios by imposing e' /e < 2.0 • 10^^, in the case ImZds < 0. 



Scenario 


A 


B 


C 


SM 


BR{Ki^ 7r°i/i>)[10-^°] 


1.3 


2.9 


11.2 


0.4 


BR{Ki^ vrOe+e") [10-11] 


2.9 


5.1 


18.2 


0.7 


BR{K+ 7r+z^P)[10-i°] 


2.0 


2.7 


4.6 


1.1 



Table 14: Upper bounds for the rare decays Ki^ — > n vu, K]^ tt e^e and tt'^vu, 
obtained in various scenarios by imposing e' /e < 3.0 • 10^^, in the case ImZ^g < 0. 



Scenario 


A 


B 


C 


SM 


BR{Kl 7r°zyi^)[10-i°] 


3.9 


6.5 


17.6 


0.4 


BR{Ki^ ^7r0e+e-)[10~"] 


7.9 


11.5 


28.0 


0.7 


BR{K+ 7r+i/z^)[10-i°] 


2.6 


3.5 


6.1 


1.1 



TT^h'i') by 1-2 orders of magnitude above the Standard Model expectations is according to 



1 174] still possible in four-generation models. 



10 Summary and Outlook 

I hope that I have convinced the students that the field of CP violation and rare decays plays 
an important role in the deeper understanding of the Standard Model and particle physics 
in general. Indeed the field of weak decays and of CP violation is one of the least understood 
sectors of the Standard Model. Even if the Standard Model is still consistent with the existing 
data for weak decay processes, the near future could change this picture dramatically through 
the advances in experiment and theory. In particular the experimental work done in the next 
ten years at BNL, CERN, CORNELL, DA^>NE, DESY, FNAL, KEK, SLAC and eventually 
LHC will certainly have considerable impact on this field. 

Let us then make a list of things we could expect in the next ten years. This list is 
certainly very biased by my own interests but could be useful anyway. Here we go: 

• The error on the CKM elements \Vcb\ and |T4b/V"cfe| could be decreased below 0.002 
and 0.01, respectively. This progress should come mainly from Cornell, i?-factories and 
new theoretical efforts. It would have considerable impact on the unitarity triangle 
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and would improve theoretical predictions for rare and CP-violating decays sensitive to 
these elements. 



The error on mt should be decreased down to ±3 GeV at Tevatron in the Main Injector 
era and to ±1 GeV at LHC. 

The measurement of non- vanishing ratio of e' je by NA31 and KTeV, excluding con- 
fidently the superweak models, has been an important achievement. The improved 
measurements of e'/e with the accuraccy of ±(1 — 2) • 10~^ from NA48, KTeV and 
KLOE should give some insight into the physics of direct CP violation inspite of large 
theoretical uncertainties. In this respect measurements of CP-violating asymmetries in 
charged B decays will also play an outstanding role. These experiments can be per- 
formed e.g. at CLEO since no time-dependences are needed. The situation concerning 
hadronic uncertainties is quite similar to e'/e. Therefore one should hope that some 
definite progress in calculating relevant hadronic matrix elements will also be made. 

More events for 'k^vv could in principle be reported from BNL already this year. 

In view of the theoretical cleanliness of this decay an observation of events at the 2-10"^'^ 
level would signal physics beyond the Standard Model. A detailed study of this very 
important decay requires, however, new experimental ideas and new efforts. The new 
efforts [113, 114| in this direction allow to hope that a measurement of Br(K~^ — > n'^i'i') 



with an accuracy of ±10% should be possible before 2005. This would have a very 
important impact on the unitarity triangle and would constitute an important test of 
the Standard Model. 

The future improved inclusive B ^s^l measurements confronted with improved 
Standard Model predictions could give the first signals of new physics. It appears that 
the errors on the input parameters could be lowered further and the theoretical error 
on Br{B Xs^) could be decreased confidently down to ±8% in the next years. The 
same accuracy in the experimental branching ratio will hopefully come from Cornell 
and later from KEK and SLAC. This may, however, be insufficient to disentangle new 
physics contributions although such an accuracy should put important constraints on 
the physics beyond the Standard Model. It would also be desirable to look for B ^d7, 
but this is clearly a much harder task. 

Similar comments apply to transitions B Xgl^l^ which appear to be even more 
sensitive to new physics contributions than B Xg^^^. An observation of — > X^/i/i 
is expected from DO and S-physics dedicated experiments at the beginning of the next 
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decade. The distributions of various kind when measured should be very useful in the 
tests of the Standard Model and its extensions. 



The theoretical status of ETl — > vr'^e+e" and of Ki^ — > fip,, should be improved to 
confront future data. Experiments at DA<I>NE should be very helpful in this respect. 
The first events of Ki^ 7r^e~^e~ should come in the first years of the next decade from 
KAMI at FNAL. The experimental status of Ki^ fifl, with the experimental error of 
±7% to be decreased soon down to ±1%, is truly impressive. 

The newly approved experiment at BNL to measure Br{Ki^ tt^vu) at the ±10% 
level before 2005 may make a decisive impact on the field of CP violation. In particular 
Ki^ — > -K^vv seems to allow the cleanest determination of ImA^. Taken together with 
TT'^ui' a very clean determination of sin 2/3 can be obtained. 

The measurement of the B^—Bg mixing and in particular of -B — > Xs^d^v and Bg^d ^jl 
will take most probably longer time but as stressed in these lectures all efforts should be 
made to measure these transitions. Considerable progress on i?^ — B^ mixing should be 
expected from HERA-B, SLAC and TEVATRON in the first years of the next decade. 
LHC-B should measure it to a high precision. With the improved calculations of ^ in 



( 4.58 ) this will have important impact on the determination of \ Vtd\ and on the unitarity 



triangle. 

Clearly future precise studies of CP violation at SLAC-B, KEK-B, HERA-B, COR- 
NELL, FNAL and LHC-B providing first direct measurements of a, /? and 7 may 
totally revolutionize our field. In particular the first signals of new physics could be 
found in the {g,fj) plane. During the recent years several, in some cases quite sophisti- 
cated and involved, strategies have been developed to extract these angles with small or 
even no hadronic uncertainties. Certainly the future will bring additional methods to 
determine a, [5 and 7. Obviously it is very desirable to have as many such strategies as 
possible available in order to overconstrain the unitarity triangle and to resolve certain 
discrete ambiguities which are a characteristic feature of these methods. 

The forbidden or strongly suppressed transitions such as — mixing and i^L — > /^e 
are also very important in this respect. Considerable progress in this area should come 
from the experiments at BNL, FNAL and KEK. 

On the theoretical side, one should hope that the non-perturbative methods will be 
considerably improved so that various Bi parameters will be calculated with sufficient 
precision. It is very important that simultaneously with advances in lattice QCD, 
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further efforts are being made in finding efficient analytical tools for calculating QCD 
effects in the long distance regime. This is, in particular very important in the field of 
non-leptonic decays, where one should not expect too much from our lattice friends in 
the coming ten years unless somebody will get a brilliant idea which will revolutionize 
lattice calculations. The accumulation of data for non-leptonic B and D decays at 
Cornell, SLAC, KEK and FNAL should teach us more about the role of non-factorizable 
contributions and in particular about the final state interactions. In this context, in 
the case of K-decays, important lessons will come from DA<I>NE which is an excellent 
machine for testing chiral perturbation theory and other non-perturbative methods. 

In any case the field of weak decays and in particular of the FCNC transitions and of CP 
violation have a great future and one should expect that they could dominate particle physics 
in the first part of the next decade. Clearly the next ten years should be very exciting in this 
field. 
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